N88- 1 1597 


5 

MASS LOSS 


Leo Goldberg 


OVERVIEW 

The first observational evidence for mass 
loss from cool stars was the discovery by 
Adams and MacCormack (1935) of blue- 
displaced components in low-excitation lines of 
metals in the spectra of certain red giant stars. 
The first systematic study of displaced lines was 
made by Deutsch (1956, 1960), who demon- 
strated that the matter forming the blue-shifted 
lines was actually escaping from the M-type 
component of the a Her system and that the 
phenomenon was common to all M giants and 
supergiants. Deutsch derived estimates of the 
mass loss for these stars, in the neighborhood 
of 10' 7 M©/yr — ten million times as great as 
the amount lost through the solar wind — and 
enough to affect the course of the late stages 
of stellar evolution in a major way. 

Weymann (1962a) carried out a detailed 
spectrophotometric study of the spectrum of 
a Ori, the first real attempt at a quantitative 
analysis of the spectrum of the outflowing mat- 
ter. He found that, despite variations of up to 
10 km/s in the radial velocity of the photo- 
sphere, the displaced lines had not varied by as 


much as 2 km/s over an interval of more than 
20 years, which implied a detached shell far 
enough away to be uncoupled from the pulsa- 
tion of the photosphere. He estimated the in- 
ner radius of the shell to be at least 14 stellar 
radii, from which the rate of mass loss was in- 
ferred to be 4.5 x 10~ 6 M©/yr. Weymann 
(1962b) also demonstrated that the mass loss 
in a Ori could not be driven by thermal gas 
pressure, as it is in the solar corona. 

After Weymann’s work and following the 
development of more sensitive detectors, in- 
terest shifted to the infrared region of the spec- 
trum. Red giants were found to emit an excess 
of infrared radiation, as compared with a 
blackbody, in a broad band of wavelengths 
centered at about 9.7 /tm. Woolf and Ney 
(1969) attributed the emission to dust grains in 
an extensive circumstellar shell and suggested 
that the grains were made of silicates contain- 
ing metallic impurities such as Fe, Mg, and A1 
(e.g., MgSi0 3 and olivine, (Mg,Fe) 2 Si0 4 ). 
From calculations of chemical equilibrium, 
Gilman (1969) showed that such molecules can 
condense in cool atmospheres of oxygen-rich 
stars (K and M stars) and identified similar 
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emissions at 11 /an as condensates of C in 
carbon-rich stars (R and N stars) and of SiC in 
stars with equal abundances of carbon and ox- 
ygen (S stars). Hoyle and Wickramasinghe 
(1962) had suggested earlier that graphite grains 
would condense in the outer regions of C-rich 
stellar atmospheres with sufficiently low tem- 
peratures. Gehrz and Woolf (1971) supposed 
that radiation pressure on the dust grains would 
propel them away from the star, dragging the 
gas along with it, and, using observed color ex- 
cesses at 1 1 /on to estimate rates of mass loss 
for more than 50 stars, they derived values in 
the range 5 x 10" 8 to 3 x 10~ 5 M© /yr. There 
was still serious doubt as to whether the mo- 
mentum acquired by the grains could be trans- 
mitted to the gas, in view of the very long mean 
free path for grain-gas collisions (Weymann, 
1962b) until Gilman (1972) showed that the gas 
and dust are momentum-coupled and not 
position-coupled (i.e., primary collisions trans- 
mit momentum from dust to gas, but secondary 
collisions between atoms (molecules) are the 
principal means by which momentum is dis- 
tributed through the gas). The process was 
shown to work for both silicate and graphite 
grains. 

The large values of mass loss inferred from 
the dust emission heightened interest in mass 
loss from red giants as an important factor in 
the late stages of stellar evolution and as a 
major source both of the enrichment of inter- 
stellar gas by heavy elements and of the inter- 
stellar dust itself. As one result, Reimers 
(1975a, 1975b) undertook a systematic program 
of analysis of circumstellar absorption lines in 
the optical spectra of more than 120 K and M 
giants and G, K, and M supergiants. In these 
and later papers (see Reimers, 1981), Reimers 
found that the stars showing circumstellar ab- 
sorption lines in their spectra, as evidence of 
high rates of mass loss, were confined to the 
upper right-hand corner of the Hertzsprung- 
Russell (HR) diagram and that terminal flow 
velocities ranged from about 5 to over 100 
km/s. The wind velocity was found to increase 
roughly as the square of the escape velocity 


from the photosphere. Except in a few cases, 
Reimers did not publish mass-loss rates derived 
for individual stars, but instead presented his 
results in the form of an interpolation formula: 

M = 4 x 1(T 13 L/gR (M©/yr) , (5-1) 

where L is the star’s luminosity, R is its radius, 
and g is the acceleration of gravity, all in solar 
units. This equation was derived from dimen- 
sional considerations, its physical significance 
being that the kinetic energy required to carry 
the escaping mass from a star to infinity is 
always the same fraction of the star’s luminos- 
ity. The constant was evaluated by the process 
of fitting the formula to 4 or 5 normal points 
derived from the empirical mass-loss rates. 
Because of uncertainties in the empirical rates, 
as well as in L, M, and R , the proportionality 
constant is generally believed to be accurate 
only within a factor of 10. Reimers’ formula 
has been the principal source of information on 
empirical rates of mass loss from red giants 
used in theoretical calculations of stellar 
evolution. 

More refined methods of analysis were ap- 
plied to a relatively small number of selected 
stars by Sanner (1976), Bernat (1977), and 
Hagen (1978). Hagen attempted to derive 
separate rates of mass loss for both the gas 
(from the optical lines) and the dust (from the 
10-fim emission feature). The absence of 
reasonable agreement among the rates derived 
by these investigators made it apparent that far 
too little is known about the structure and com- 
position of red giant atmospheres, as will be- 
come evident further on. 

Observations with the International Ultra- 
violet Explorer (IUE) and with the High Energy 
Astrophysical Observatory (HEAO) have illum- 
inated further the nature of the division in the 
HR diagram between stars with high and low 
rates of mass loss. The short-wavelength limit 
of observation with IUE is about 1 170 A, which 
reaches lines formed in the chromosphere/co- 
rona transition region, but not the coronal lines 
themselves. What at first appeared to be a 


246 


sharp line dividing stars with hot coronas from 
those with cool massive winds (Linsky and 
Haisch, 1979) is now seen as a gradual transi- 
tion, with some stars in the boundary zone 
showing both types of spectra (Hartmann et al., 
1980; see also Dupree, 1980; Reimers, 1981). 
Stars near the boundary tend to show temporal 
variations in the mass flows. 

New techniques for observing at infrared, 
millimeter, and centimeter wavelengths are cur- 
rently leading to dramatic advances in knowl- 
edge of the rates and evolutionary implications 
of mass loss from highly evolved stars. Obser- 
vations with very high resolution of the vibra- 
tion/rotation line spectra of CO, including both 
fundamental and overtone bands, have been ef- 
fective in probing the structures of mass flows 
from red giants, especially Miras (Hall, 1980). 
Emission from circumstellar shells in the 1-0 
and 2-1 transitions of CO at millimeter wave- 
lengths can give direct information on stellar 
velocities, shell expansion velocities, shell sizes, 
and, with appropriate modeling, estimates of 
mass loss (Knapp et al., 1982; Knapp and Mor- 
ris, 1984). Flux densities of radio continuum 
radiation measured with the very large array 
(VLA) may now be used to infer the tempera- 
ture and extent of chromospheres in red giants. 
The VLA has also been used to map the a Sco 
system (Hjellming and Newell, 1983), in which 
the creation of an H II region in the expanding 
shell of the M star by the ultraviolet (UV) radia- 
tion of the B star permits a determination of 
the rate of mass loss from the M star. Finally, 
maser radiation emitted by H 2 0, SiO, and OH 
has been shown in principle to yield informa- 
tion on the structure of various regions of late- 
type atmospheres and on the enormous extent 
of their circumstellar shells (see Elitzur, 1981). 
The so-called OH/IR stars offer the possibility 
of studying very luminous, long-period Mira 
variables which are not seen in visible light. 

The observational study of mass loss from 
red giants in globular clusters was initiated by 
Cohen (1976), who interpreted emission asso- 
ciated with the Ha absorption line as arising 
by recombination in an expanding ionized shell 
of hydrogen. The same method has been ex- 


tended to large numbers of stars by Mallia and 
Pagel (1978), Peterson (1981, 1982), and espe- 
cially by Cacciari and Freeman (1981, 1983), 
all of whom derived rates in the neighborhood 
of 10' 8 M 0 /yr. Although rates of this order 
of magnitude seem to satisfy the requirements 
of stellar evolution, there is a fundamental dif- 
ficulty in the assumption that the presence of 
Ha emission is by itself evidence of mass loss 
(Reimers, 1981; Dupree, et al, 1984). 

In addition to thermal gas pressure and ra- 
diation pressure on dust grains, a number of 
other interesting mechanisms have been pro- 
posed for both steady and episodic mass loss 
(e.g., wave pressure, shock waves, magnetic 
reconnection, dynamical and pulsational insta- 
bility, and pulsational mode switching in 
Miras). The probability is high that the dust 
mechanism operates in the coolest late-type 
stars, but it is not at all clear which mechanisms 
are most relevant in other stars. 

Mass loss from red giants plays a crucial role 
in several aspects of stellar evolution (see Ren- 
zini, 1981a; Iben, 1981b). Mass loss from stars 
of intermediate mass between 1 and 8 M© can 
make an important contribution to the enrich- 
ment of the interstellar medium. Moreover, in 
the absence of mass loss, all stars with masses 
exceeding 1.4 M© would become supernovae, 
but if enough mass is lost while the star is on 
the asymptotic giant branch (AGB) of the HR 
diagram, intermediate-mass stars may evolve 
into planetary nebulae and white dwarfs. Par- 
ticularly exciting is the recognition of long- 
period OH/IR stars as possible progenitors of 
planetary nebulae. Mass loss from interme- 
diate-mass stars may also affect the formation 
of interstellar dust grains, Mira variables, and 
carbon stars. The evolution of low-mass stars 
(mass less than 1.4 M© ) has been extensively 
studied in globular clusters. In the absence of 
mass loss, one finds discrepancies between the 
observed and calculated luminosity functions 
of stars on the horizontal branch (HB) and 
AGB, which can be removed by postulating 
mass losses of about 0.2 M© on the red-giant 
branch (RGB) and about 0. 1 A/© on the AGB. 


247 


EMPIRICAL RATE 
DETERMINATIONS 

UV-Optical Spectra 0.2-5 /im 

Single Stars. Figure 5-1 shows small regions of 
the spectra of several red giants and supergiants 
centered on the resonance lines of Mg + , Ca + , 
Na, Mn, and K. The profiles are more or less 
typical of the resonance lines of neutral and 
singly ionized metals in the visible and near-UV 
spectral regions of cool stars showing mass loss. 
Figure 5-2 shows a tracing of the R1 line of the 
1-0 band of CO near 4.65 jim, which is similar 
in appearance to the potassium line. The cores 
of the Mg II and Ca II lines are strongly satu- 
rated and contaminated by chromospheric 
emission, but the other lines show the P Cygni 
profiles characteristic of resonance scattering 
in an extended expanding atmosphere, although 
distorted by the underlying, broad photospheric 
line. The absorption components are formed 
along the line of sight to the star and are dis- 
placed to the blue by the amount of the expan- 
sion velocity, whereas the emission is undis- 
placed relative to the center-of-mass velocity. 
Notice that in a Ori, the CO and K I cir- 
cumstellar (CS) lines are double, with an ap- 
parent separation of about 6.5 km/s, implying 
the existence of two discrete shells. Curiously, 
in this star, the doubling is also seen in the D- 
lines of Na I, but not in any other lines. The 
absence of shell components among lines with 
excitation potentials greater than a few tenths 
of a volt has been interpreted as implying shell 
temperatures less than 1000 K (Weymann, 
1962a). 

Spectral line profiles are used for the deriva- 
tion of mass-loss rates with the aid of the equa- 
tion of continuity: 

M = 4t • r 2 * n(r) • • v(r) , (5-2) 

where r is the radius of a volume element with 
hydrogen number density n(r) and outward 


velocity u(r), fi is the mean molecular weight, 
and m H is the mass of the H atom. This equa- 
tion implies steady mass loss with spherical 
symmetry. In most cases, the radial dependence 
of velocity and density cannot be determined 
without additional assumptions. Since the cir- 
cumstellar absorption lines appear to be 
relatively sharp, the velocity along the line of 
sight is usually assumed to be constant, which 
implies an inverse square law for the density. 
Moreover, the formation of P Cygni profiles 
signifies that most of the emission is occurring 
far from the star, and the radial velocities of 
the displaced components are constant even in 
stars with variable photospheric velocities. 
Therefore, the circumstellar shells are assumed 
to be detached from the stellar photospheres 
and to be far enough away to be insensible to 
their motions. With these two assumptions, 
Equation (5-2) becomes: 

M = 4tt * N- ii- m H - R o * V , (5-3) 

where TV is the column density of hydrogen gas 
along the line of sight, R is the inner radius 
of the detached shell, and V is the expansion 
velocity, assumed to be constant. The deriva- 
tion of M from Equation (5-3) proceeds by 
steps as follows. 

First, the P Cygni profile must be detached 
from the background photospheric line, which 
can be done with varying degrees of difficulty. 
Two quantities are extracted from the line pro- 
file: the column density of atoms or molecules 
in the ground level of a particular transition and 
the expansion velocity. The earlier work 
(Weymann, 1962a; Reimers, 1975a) assumed 
line formation in a static, plane-parallel at- 
mosphere, but more recent work has made use 
of relatively exact calculations by Kunasz and 
Hummer (1974) of radiative transfer in an ex- 
panding spherical atmosphere and of relatively 
sophisticated non-LTE calculations of ioniza- 
tion equilibrium (Sanner, 1976; Bernat, 1977; 
Hagen, 1978). On the average, column densities 
derived by Bernat (1977) for a Ori are smaller 
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Figure 5 - 7 . Circumstellar resonance line profiles: (a) K I in Ori; (b) Mg II in a Ori (Bernat 1976); 
(c) Na 1 5890 A in a Her (left); Ca II K in 0 Peg (Hagen, 1980) (right); (d) Mn I in a Ori . The vertical 
line labeled “em” in (a) marks the center of the emission component of the P Cygni profile. 
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Figure 5-2. The R1 line of the 1-0 band of CO 
in the spectrum of a Ori . The two shell com- 
ponents , s ] and s r are superposed on the 
broad photospheric line labeled p. The strong 
telluric component is shown to the right (Ber- 
nat etai, 1979). The vertical line labeled “em” 
marks the center of the emission component of 
the P Cygni profile. 

than those of Weymann (1962a) by a factor of 
3 to 4. The expansion velocity, V 9 is obtained 
by direct measurement of the displacement of 
the absorption line relative to neighboring 
photospheric lines. The emission on the red- 
ward side of the absorption line causes a 
spurious additional shift of the line, amounting 
to an error of about 25 percent (Bernat, 1976). 
If the velocity of the photosphere is variable, 
the expansion velocity is the shell velocity minus 
the center-of-mass velocity. 

The column density of the given atom in a 
given stage of ionization must then be con- 
verted into, first, the column density of the ele- 
ment in all stages of ionization and, second, the 
column density of hydrogen. Most of the res- 
onance lines seen in the visible spectrum, with 
the exception of those of Ca II, Sr II, and Ba 
II, belong to neutral metals such as Fe, Cr, Ni, 
Mn, and K. Since the metals are predominant- 
ly singly ionized, a theory of ionization is need- 
ed to derive the abundance of the ion. This step 
is omitted for Sr II and Ba II. The total col- 
umn density then follows from the assumption 


that the abundances are solar-like, which may 
be reasonable for most of the metals but is 
questionable for Sr and Ba. The assumptions 
made for CO in O-rich stars are that all car- 
bon atoms are locked up in CO, which seems 
reasonable, and that the C/H ratio is the same 
as in the Sun, which is a pure guess. 

Applications of the theory of ionization 
have led to severe contradictions in the results. 
As a test of the theory, Hagen (1978) derived 
column densities of Sr II and Ba II in eight stars 
and then used solar abundances and ionization 
theory to predict column densities for several 
neutral metals. Figure 5-3 shows the logarithm 
of the ratio of the observed column density 
N(obs) to the predicted column density N(pred) 
plotted against the element in order of decreas- 
ing ionization potential. The discrepancy of 
four orders of magnitude is hard to explain. 
Part of it apparently results from the steep drop 
of the assumed blackbody radiation field below 
2000 A, and part may be caused by the deple- 
tion in grains of elements with higher ioniza- 
tion potentials. Figure 5-3 suggests that neutral 
elements of low ionization potential are much 
more abundant than the theory predicts and 
raises a question as to whether Ba and Sr are 
indeed predominantly singly ionized, at least in 
the cooler M giants. 

Hagen (1978) used only the lines of Sc II, 
Sr II, and Ba II to derive column densities of 
hydrogen in eight red giants. Bernat (1977) 
derived hydrogen column densities in the four 
red supergiants— a Ori, a Sco, a Her, and p 
Cep— by averaging the results obtained from 
8 to 9 metallic atoms, both neutral and singly 
ionized. Although the scatter is large, ap- 
proaching a factor of 50, there is no systematic 
difference between hydrogen column densities 
derived from neutral and ionized species, 
respectively. On the other hand, the column 
densities inferred from the Ti II lines at 3300 
A are too small relative to those from Sr II and 
Ba II by factors of 7 to 20, if solar relative 
abundances are used and if the singly ionized 
stage is assumed to be predominant. 
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Figure 5-3 . The logarithm of the ratio of the 
observed column density of neutral atoms to 
the column density predicted through the use 
of ionization theory and solar abundances plot- 
ted against element in order of decreasing 
ionization potential. The numbers refer to dif- 
ferent stars as follows: 1 = X Her , 2 = RX 
Boo , 3 — HD 207076 , 4 = R Leo , 5 ~ a Sco, 
6 = a Ori , 7 - a Her , 8 = fi Cep (Hagen, 
1978). 

A very large error in M may arise because 
there is often no satisfactory method available 
for the determination of R , the inner radius 
of the shell; different investigators have used 
values differing by one or two orders of 
magnitude for the same star. Weymann (1962a) 
proposed a method for deriving the inner radius 
based on the idea that the cores of the infrared 
triplet lines of Ca II, 8542 and 8662 A, might 
show small displaced absorption components 
formed in the expanding shell. Since the lower 
levels of these lines should be populated by 
radiative excitation, the equivalent width of the 
feature ought to be inversely related to the shell 
radius. The method has been applied to several 


stars by Weymann (1962a), Sanner (1976), and 
Bernat (1977), and although it appears to work 
in some stars, it leads to widely discrepant 
results for a Ori, for which the lines are so deep 
that the signal to noise tends to be relatively 
small at the very line center (Goldberg, 1979). 
Another method of estimating R q is based on 
model calculations of off-limb profiles of K I 
emission (Bernat, 1976). When the line of sight 
is interior to the inner radius, the calculated 
profiles are asymmetric, and the observation of 
such asymmetries in Betelgeuse implies an in- 
ner radius of about 100 stellar radii, but the 
calculations are strongly model-dependent. 

There is in fact some question as to whether 
the shell has a well-defined inner radius, as is 
usually assumed in mass-loss determinations. 
No matter what the mechanism for initiating 
and maintaining the flow may be, there must 
be a region in which the gas accelerates to its 
terminal velocity and in which the density is 
relatively high. Such a region would be difficult 
to identify from a P Cygni profile because, al- 
though the absorption component might be 
strongly asymmetrical, the asymmetry would be 
camouflaged by the redward-shifted emission 
component. The superposition of a P Cygni 
profile on a photospheric line would make such 
an asymmetry doubly difficult to detect. On the 
other hand, subordinate lines would not be ex- 
pected to display strong redward emission, since 
photons absorbed in such lines tend to be re- 
emitted in resonance lines. Such asymmetric 
subordinate lines are indeed observed in the 
spectrum of a Ori, the most prominent among 
them being the Fe I lines at 4140 A and the Ca 
II IR triplet (Goldberg, 1979), both of which 
have excitation potentials of about 1.5 eV. The 
8498 and 8542 A lines of Ca II are shown in 
Figure 5-4, where they are plotted on a helio- 
centric velocity scale. Notice that the line cores 
are highly asymmetric in the sense that the 
radial velocity decreases from about 21 to 22 
km/s, which is the systemic velocity of the star, 
in the upper part of the core, to about 13.5 
km/s at the bottom of the 8542 A line. More- 
over, over a period of several years, during 
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Figure 5-4. Profiles of the 8498 and 8542 A lines 
of Ca II in a Ori made with the Fourier 
transform spectrometer of the McMath Tele- 
scope at Kitt Peak National Observatory, The 
lines are plotted on a velocity scale. The line 
bisectors are indicated by dots (8498 A) and 
crosses (8542 A), with the corresponding 
heliocentric velocities shown to the left. 


which the mean velocity of the photosphere 
decreased by about 5 km/s, the velocities in 
these profiles remained stationary (Goldberg, 
1979). The terminal velocity shown by the ab- 
sorption components of the resonance lines is 
about 1 1 km/s, but if a correction is made for 
the spurious shift caused by their companion 
emission lines, the velocity would agree with 
that measured at the center of the 8542 A line. 
Thus, the line is formed in. an atmospheric 
region which is far enough away to be un- 
coupled from the motion of the photosphere, 
which is an argument in favor of the detached 
shell model. It is probable that the absorption 
components of the resonance lines would show 
the same behavior if they could be separated 
from their emission components. 


The first high-resolution observations of cir- 
cumstellar CO were made in the spectrum of 
a Ori by Bernat et al. (1979) with the Fourier 
transform spectrometer (FTS) at the Kitt Peak 
National Observatory (KPNO) 4-meter tele- 
scope. Numerous lines of the fundamental 
vibration/rotation band at 4.6 gm were ob- 
served, including several lines of the isotopic 
molecule, 13 CO. Subsequently, Bernat (1981) 
observed the band in nine other red giants. 
Figure 5-2 shows a typical CO line in the spec- 
trum of Betelgeuse, which consists of four com- 
ponents: a telluric line, a broad photospheric 
line, and two narrow circumstellar lines with 
heliocentric velocities of 11.0 and 4.8 km/s, 
respectively. The same velocity components are 
observed in K I and probably also in the D-lines 
of Na (Goldberg et al., 1975), but in no other 
lines. The relative abundances 12 C/ 13 C were 
found to be about 6 for both components, 
which suggests two shells with a common origin 
in highly processed matter from the stellar in- 
terior and thus rules out an interstellar origin 
for the second component. After showing that 
collisional excitation dominates over radiative 
excitation, Bernat et al. (1979) identify the ex- 
citation temperature with the kinetic 
temperature, assume that the gas and dust are 
in thermal equilibrium, and derive distances of 
150 and 2000 stellar radii, respectively, for the 
two shells. The further assumptions of cosmic 
composition and the complete association of 
carbon into CO lead to mass-loss rates of 5 x 
10" 6 and 3.2 x 10' 6 M© /yr for the near and 
far shells, respectively. 

The results of Bernat’s (1981) survey of 
other red giants for circumstellar CO are shown 
in Table 5-1. In addition to these stars, three 
other stars— a Boo, a Tau, and a Sco— showed 
no detectable CO, whereas the CO shell lines 
in RX Boo were found to be too strongly sat- 
urated for analysis, a Boo and a Tau are too 
early in spectral type to show circumstellar 
lines, and as Bernat points out, any CO in the 
shell of a Sco would be dissociated by the com- 
panion B star. Table 5-1 shows that the CO ab- 
sorption features are multiple in many stars, 
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Table 5-1 

Observations of CO in M Giants and Supergiants (Bernat, 1981) 


Star 

Spectral 

Type 

Component 

(km/s) 

T 

exc 

(K) 

N (CO)/v (Dop) 
(cm -2 km -1 s) 

N (H) 
(cm' 2 ) 

N (Dust) 3 
(cm' 2 ) 

119 Tau 

M2 lab-lb 

-9 

200 ±150 

7.5 + 1 5 

3.4 + 21 b 

— 

fi Cep 

M2 la 

-8 

100±10 

7.5+15 

7.3 + 21° 

— 



-13 

270 ±60 

2.0+16 

— 

— 



-19 

1 00 ± 1 5 

1.8+16 

— 

4.0 + 21 d 



-38 

60 ± 4 

1.1 + 17 

— 

— 



-47 

1 00 ± 40 

3.0+15 

— 

— 

(3 Peg 

M2— 311— III 

-6 

90 ±30 

1.3+16 

2.7 + 20 b 

— 

q Per 

M4 ll-lll 

-2 

90 ±20 

6.5+15 

2.7 + 20 b 

— 

a Her 

M5 Ib-ll-f 

-13 

2 50 ± 60 

4.2 + 16 

6.6 + 21 d 

2.0 + 20 d 



-25 

550± 670 

6.7+15 

— 

— 

SW Vir 

M7 III 

-6 

1 30 ± 20 

2.9+16 

— 

— 



-9 

1 30 ± 1 5 

2.3+16 

— 

— 

X Her 

M6e 

-8 

1 1 0 ± 20 

4.3+16 

2.5 + 20 d 

4.0 + 21 d 

W Hya 

M7-9e 

-5 

300 ±90 

4.8+16 

2.5 + 20 d 

6.0 + 21 d 



-13 

1 20 ± 20 

3.7+15 

— 

— 


a 

b 


c 

d 


N(Dust) is the H column density derived assuming cosmic abundances and complete condensation of silicates. 
Sanner (1976). 

Bernat (1977). 

Hagen (1978). 


which Bernat now attributes to episodic mass 
loss rather than to the formation of discrete 
clumps by instabilities in a steady flow as was 
postulated for Betelgeuse (Bernat et al., 1979). 
Column 4 lists the excitation temperatures, in- 
terpreted as kinetic temperatures. The Doppler 
velocity could not be independently determined 
and therefore only the ratio of CO column den- 
sity to Doppler velocity can be inferred, as given 
in column 5. Column 6 contains the hydrogen 
column densities derived in previous investiga- 
tions of optical atomic spectra, and column 7 
lists the hydrogen column density derived from 
the \0-fim silicate feature on the assumption 
that all Si atoms are condensed on grains and 
that the Si/H ratio has the solar value. The 
absence of a correlation between CO column 
densities and either of the two hydrogen column 
densities is, in Bernat’s view, a major argument 
in favor of episodic mass loss, since the separate 
shells would be formed under different physical 
conditions. 


Bernat’s final conclusion is that meaningful 
mass-loss rates cannot be determined from ob- 
servations of nonsteady mass loss unless aver- 
ages are taken over several episodes. We shall 
come back later to the topic of episodic mass 
loss and comment here only that the absence 
of correlations between the various column 
densities is not unexpected in view of: (1) the 
large scatter of the determinations from differ- 
ent ions in the same star, and (2) the many 
sources of error in the derivation of column 
densities from the 10-/un emission feature. 

Binary Stars. The errors introduced by un- 
certain knowledge of the degree of ionization 
and of the inner shell radius may be avoided 
in principle by observing circumstellar lines of 
elements in a predominant stage of ionization 
(e.g., Fe II or Ti II) in the spectrum of a hot 
visual companion of a red giant (Reimers, 
1981). If the geometry of the binary system is 
known, no assumption need be made about an 
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inner shell radius. Stars studied in this way from 
the ground, from a balloon, and from IUE are 
a 1 Her (Reimers, 1977a) and a Sco A (Kudrit- 
zki and Reimers, 1978; van der Hucht et al., 
1980; Bernat, 1982). Reimers (1977a) derived 
hydrogen column densities from six shell lines 
of Sr II (4077 A) and Ti II (3300 A) in the spec- 
trum of a 2 Her (G5 III), the visual companion 
of a ] Her (M5 II), and obtained a mass-loss 
rate of 1.1 x 10" 7 M 0 /yr, which is about 
one-sixth of Bernat’s (1977) value derived from 
lines of neutral metals at visible wavelengths in 
the spectrum of a 1 Her. 

The binary-star approach has been exten- 
sively applied to the a Sco system. Kudritzki 
and Reimers (1978) derived a rate of 7 x 10“ 7 
Mq /yr for mass loss from a Sco A (Ml. 5 lab), 
using the H and K lines of Ca II and the Ti II 
lines in the spectrum of a Sco B (B2.5 II). Van 
der Hucht et al. (1980) observed a variety of 
lines between 2000 and 3000 A in the spectrum 
of a Sco B with the balloon-borne ultraviolet 
spectrometer (BUSS), giving greatest weight to 
the Zn II and Cr II lines, and obtained a rate 
equal to 7.1 x 10~ 6 Mq /yr. It is possible that 
the Zn II lines are primarily interstellar 
(Reimers, 1981). Finally, Bernat (1982) has 
made the most thorough investigation of mass 
loss from a Sco A, with the IUE satellite, by 
observing 90 lines of about a dozen ions in the 
wavelength range 1190 to 2631 A. The result 
is that the column densities of hydrogen, de- 
rived from 10 different ions, predominantly 
singly ionized, show a spread of a factor of 30. 
No single explanation— observational error, 
shell chemistry, interstellar contributions— 
appears to be satisfactory. An unweighted 
mean over all ions gives a rate of 6.4 x 10~ 6 
M 0 /yr. 

In principle, the binary-star approach may 
be extended to a much larger number of stars 
by IUE observations of spectroscopic binaries 
consisting of a red-giant star and a B-type com- 
panion. Observations in the ultraviolet effec- 
tively separate the two stars. The first attempts 
along these lines were made for 32 Cyg (Stencel 
et al., 1979) and f Aur (Chapman, 1981), but 


neglected to take into account the nonspherical 
nature of the line-transfer problem (Reimers, 
1981). The angular separation of such objects 
is so small that radiation from both the B star 
and the circumstellar shell pass through the 
spectrograph slit, resulting in P Cygni profiles 
that must be analyzed in the context of a non- 
spherical three-dimensional line-transfer prob- 
lem. The necessary theoretical framework has 
been developed by Hempe (1982) and applied 
to the analysis of a number of objects, in- 
cluding f Aur, 32 Cyg, and 31 Cyg (Che et al., 
1983), <5 Sge (Reimers and Schroder, 1983), and 
Boss 1985 (Che and Reimers, 1983). For the K 
supergiant, f Aur, Che et al. (1983) derived M 
= 0.63 x 10~ 8 Mq/ yr and v = 40 km/s, 
as compared with Chapman’s (1981) values of 
M - 2 x 10" 8 Mq /yr and v - 100 km/s, 
respectively. For 32 Cyg, Che et al. found 
Ki = 2 x 10~ 8 Mq /yr and v =60 km/s, 
as compared with Stencel et al.’s (1979) 4 x 
10“ 7 A/ 0 /yr and 65 km/s, respectively. 

Mira Stars. Optical studies of mass loss from 
Mira stars have benefited greatly from obser- 
vations of the infrared vibration/rotation bands 
of CO. The star x Cyg is a good example, in 
which a series of high-resolution spectra, taken 
with the FTS at the KPNO 4-meter telescope 
may have provided new insights into the mass- 
loss mechanism for red-giant stars (Hinkle et 
al., 1982). The spectra covered more than three 
cycles of the visible light variation and included 
the first and second overtone bands as well as 
the fundamental. Analysis has revealed four 
sets of lines, separated according to their veloc- 
ity characteristics, and apparently arising from 
four different atmospheric regions: (1) a pul- 
sating photospheric region through which a 
shock front passes near maximum light; (2) a 
cool gas (T - 300 K) expanding at 7.5 km/s 
relative to the center of mass, which is identified 
with the circumstellar shell; (3) a hotter region 
(T - 800 K), stationary with respect to the 
center of mass; and (4) a still hotter gas (T - 
1500 K) falling inward at 8 km/s. The column 
density of the stationary layer varies on a time 
scale larger than the optical period of 410 days. 
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For example, the lines were absent in 1975, be- 
came quite pronounced 7 months later, and 
then weakened steadily over the next 3 years. 
It is conjectured that this layer contains gas 
ejected during particularly violent oscillations 
and is somehow supported by dissipated shock 
energy. If the stationary layer is heated by stel- 
lar photospheric radiation, its 800 K tempera- 
ture implies a distance of 10 /?*, or 1 .7 x 10 14 
cm, where dust grains can form. Roughly the 
same value, 10 14 cm, was measured by in- 
frared speckle imaging techniques (Mariotti et 
al., 1983). Thus, it can provide a reservoir both 
for the dust-driven mass loss and the infalling 
gas. A similar layer appears to be present in the 
10 Miras which have been observed thus far by 
Hinkle et al. (1982), who suggest that it may 
also be present in most late-type stars. 

Red Giants in Globular Clusters. It has been 
known for some time (for a summary, see Ren- 
zini, 1981a) that a mass loss of about 0.2 M© 
during the RGB (red-giant branch) phase, and 
an additional 0.1 A/© during the AGB (asymp- 
totic giant branch) phase are needed to bring 
the calculated distribution of stars on the HR 
diagrams of globular clusters into agreement 
with observations. In the absence of experimen- 
tal evidence for mass loss from stars in globular 
clusters, theorists incorporated into their cal- 
culations (Fusi-Pecci and Renzini, 1976; Ren- 
zini, 1977) Renners’ (1975a, 1975b) parametric 
formula for mass loss, Equation (5-2), modified 
by an adjustable parameter, 77 , which is sup- 
posed to reflect the uncertainty in the empirical 
rates: 

M = -4 X 10 “ 13 X rj • L/gR. (5-4) 

This relationship was derived for Population 
I stars, and the scatter around it suggests that 
r] lies somewhere between 0.3 and 3.0. It is 
therefore remarkable that, within rather close 
limits, the Reimers formula with 77 = 0.40 
clears up both the HB and AGB discrepancies, 
despite the fact that it is based on empirical 


rates for relatively few Population I stars, 
which are themselves uncertain by a factor of 
10 . 

Observational evidence for mass loss from 
red giants in globular clusters was first offered 
by Cohen (1976) in the form of emission assoc- 
iated with the Ha absorption lines in four red 
giants belonging to three globular clusters. 
Although emission has been observed on the 
blue wing in some cases, it has also been seen 
only on the red wing in other stars and on both 
wings in still others. Cohen’s observations have 
been extended to a large number of stars in 
more than a dozen clusters by Mallia and Pagel 
(1978), by Peterson (1981, 1982), and espe- 
cially by Cacciari and Freeman (1983). Cohen 
(1976) postulated that the emission arises from 
a circumstellar shell as a result of mass loss. As 
a model for estimating mass loss, Cohen took 
a completely ionized shell, expanding at a con- 
stant rate of 45 km/s from an inner shell radius 
R s = 2 R* to infinity. The envelope was 
assumed to be optically thick in the Lyman lines 
and optically thin in the Balmer lines. With 
these assumptions, the total Ha emission from 
the star is written as: 

00 

E(Ha) = 47 t • hv • • a 32 R* j dr/r 2 , (5-5) 

Jr 

S 

in ergs/s, where n Q is the density of hydrogen 
at R $y and a 32 is the effective recombination 
coefficient for a temperature of 10000 K. 

The rate of mass loss is then 

M=4 »•/?•>»„• (5-6) 

where n o is to be substituted from Equation 
(5-5), and m R is the mass of the hydrogen 
atom. Numerically, with the equivalent width 
W substituted for F(Ha), the rate of mass loss 
becomes (Mallia and Pagel, 1978): 

rif = 2.76 X 10 " 4 V 

exp 

(5-7) 

. /?* (R s • W) v2 • exp(-l.l/T 4 ) . 
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The radii are expressed in solar radii, V gxp 
is in km/s, W is in A, T is in units of 10 4 K, 
and A^f is in Mq /yr. Note that, when expressed 
in terms of the observed quantity, W, the 
equivalent width, rif, is proportional to R* 2 , 
and not to RY 2 as Cohen (1976) supposed. 

All investigators have followed Cohen’s 
(1976) methodology in calculating rates of mass 
loss. Cohen derived a minimum rate of 2 x 
10~ 9 Mq /yr, which, when combined with a 
red-giant lifetime of 1 x 10 8 years, gave the 
desired loss of 0.2 M © on the RGB. Mallia and 
Pagel (1978) found much larger rates, up to 6 
x 10' 8 Mq /yr, principally because they used 
larger stellar radii, but they also believed that 
their results were compatible with the re- 
quirements of HB morphology. In her study of 
seven luminous giants in three clusters, using 
high spectral resolution, Peterson (1981) found 
that the Na D and Ha lines in some stars had 
blue-shifted asymmetric cores similar to those 
observed in M field giants and supergiants 
(Boesgaard and Hagen, 1979). The shifts 
amounted to about 6 km/ s in Ha and 1 1 km/ s 
in Na D. Peterson (1981) notes that emission 
is not usually seen in field giants with asym- 
metric cores. Boesgaard and Hagen (1979) find 
that 50 percent of stars with spectral classes MO 
III to M7 III show asymmetric Ha profiles with 
core blue shifts of a few km/s. 

In a recent extensive survey of 143 red giants 
in 12 globular clusters, Cacciari and Freeman 
(1983) detected Ha emission in about one-third 
of the stars brighter than log L = 2.9, or M bgl 
< -2.5, with clear evidence of variability. More- 
over, the apparent narrowing of Ha absorption 
in fainter stars suggests the presence of weaker 
emission down to log L = 2.3, or M bol = -1. 
Mass-loss rates calculated from Equation (5-7) 
are between 0.7 and 6.3 in units of 10 8 
M q /yr. The required total mass loss on the 
RGB is attained if it is assumed that the average 
rate of loss from M bol = -1 to the tip of the 
RGB is 10' 8 Mq /yr over an evolutionary time 
of 2 x 10~ 7 years (Sweigart and Gross, 1978). 

Despite these estimates, however, the occur- 
rence of mass loss in globular cluster stars on 


the RGB in an amount sufficient to satisfy 
evolutionary requirements is still unproven 
(Reimers, 1981; Dupree et al„ 1984). Neither 
Ha emission nor the dependence of the emis- 
sion on luminosity is evidence of mass loss. 
Model calculations of Ha line profiles in metal- 
deficient stars by Dupree et al. (1984) show that 
emission wings can arise in warm chromo- 
spheres with mass-loss rates 2 x 10' 9 
Mq /yr. Moreover, chromospheric activity 
might well increase as stars ascend the red-giant 
branch. The assumption that the emission is 
formed in a fully ionized gas by recombination 
is purely arbitrary, and uncertainties in the 
assumed stellar radii represent another source 
of error, which has already led to a disagree- 
ment of an order of magnitude (Cohen, 1976; 
Mallia and Pagel, 1978). The observation of 
asymmetric blue-shifted line cores (Peterson, 
1981, 1982) implies that the outer envelopes of 
some red giants in globular clusters are ex- 
panding, but not necessarily that matter is 
escaping. 

Infrared Dust Emission 

Emission by silicates in a broad band 
centered at 9.7 ^m is a common occurrence in 
the spectra of late-type O-rich giants and 
supergiants. In the carbon stars, a similar but 
narrower band near 1 1 /un is probably radiated 
by SiC, while a third species of unknown com- 
position contributes to the continua of some, 
but not all, C and S stars. The silicate feature 
is the most useful indicator of mass loss in the 
middle infrared. In the far infrared, thermal 
continuum emission near 400 fim by silicate and 
carbon grains is the most promising determi- 
nant of dust mass loss from highly evolved 
stars. We consider first the silicate emission at 
10 fim and next the continuous emission at sub- 
millimeter wavelengths. 

Emission by Silicates at 10 Microns. Figure 5-5 
shows the excess emission observed in the 
10-/iin region of red-giant spectra, as compared 
with blackbody emission at the stellar 
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temperature. Basically, one tries to extract the 
optical depth of the dust grains at 10 fim from 
the amplitude of the excess emission, bearing 
in mind that the amount radiated also depends 
on the temperature of the dust. The dust 
temperature is determined by the composition 
of the grains and their absorbing properties at 
the shorter infrared wavelengths where most of 
the stellar radiation is concentrated. The radial 
optical depth T(10/un) is then used to calculate 
the column density of dust, 

7(10 pm) = a 2 • Q(10 pm) • N e , (5-8) 

where a is the radius of a grain, Q is the ab- 
sorption efficiency at 10 pm, and N is the col- 
umn density of silicate grains per cm 2 . If 



Figure 5-5. Spectral energy curves for K and M 
supergiants in the 2- to 3-pm region. The ex- 
cess emission near 10 pm is clearly seen. For 
comparison, curves for a 3300 K blackbody and 
for free-free emission are drawn on a relative 
scale (Merrill, 1977). 


values of a and O are known or assumed, N 

C 

may be found from Equation (5-8). Next, the 
column density of hydrogen may be inferred 
from: (a) the assumed grain composition, (b) 
a judgment as to the fraction of a given element 
that may have condensed on the grains, and (c) 
an estimate of the abundance of that element 
relative to hydrogen. Finally, the mass-loss rate 
is calculated from Equation (5-3) with the aid 
of expansion velocities derived from optical 
spectra and estimates of the inner shell radius 
by one means or another. 

The most uncertain step in the analysis is the 
derivation of r(10 fim) from the observed IR 
emission, which requires knowledge of half a 
dozen other parameters, and is therefore 
strongly model-dependent (Hagen, 1978): 

1. The dust temperature distribution is 
determined by the equilibrium between 
the absorption of stellar radiation at 
wavelengths smaller than 5 ^m and the 
emission at longer wavelengths. 
Knowledge of grain compositions and 
absorption efficiencies are of critical im- 
portance. In discussions of grain 
temperatures, a distinction is usually 
made between so-called 4 ‘clean ' 9 silicates, 
which are relatively pure and have highly 
developed crystalline structures, and 
“dirty” silicates, which have disordered 
lattices (e.g., natural hydrous silicates in 
carbonaceous chrondrites) or are coated 
with impurities such as graphite. Clean 
silicates are transparent to short infrared 
wavelengths, whereas the dirty silicates 
are strong absorbers at these wave- 
lengths. Jones and Merrill (1976) showed 
that pure silicate grains reach an equi- 
librium temperature that is too low to 
reproduce the observed 10-/an feature. 
The introduction of impurity (dirty 
silicates) raises the short-wave efficien- 
cy and warms the grains. 

2. The absorption efficiencies are well 
known for relatively few minerals, and 
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although they are reasonably indepen- 
dent of composition at 10 jim, they dif- 
fer from one mineral to another by 
several orders of magnitude at wave- 
lengths less than 5 jim, where the stars 
radiate most of their energy. The tem- 
perature distributions are correspond- 
ingly uncertain. 

3 . For given optical depth, the size of the 
shell, particularly the inner radius, deter- 
mines the flux at 10 jim. 

4 . The density distribution for dust grains 
is another important input. It is usually 
assumed that the density falls off as r" a , 
where a = 2 for constant velocity. Hagen 
(1978) prefers a = 1.5, but nevertheless 
assumes constant velocity. 

5 . The small-grain approximation (Shawl, 
1975) is usually chosen to fix the grain 
size, with a = 0.05 to 0.1 jim. In this 
range, a « X, and Q(abs) oc a. The dust 
temperature is independent of grain size, 
since both the energy absorbed per unit 
mass and the energy emitted are propor- 
tional to the grain volume. 

6. Hagen (1978) took partial account of 
radiative-transfer effects by including a 
correction for self-absorption of starlight 
by the dust, but a more exact treatment 
is required when r(10 jim) is greater than 
unity (Rowan-Robinson and Harris, 
1983a). 

The extensive survey of IR circumstellar 
emission by Gehrz and Woolf (1971) was based 
on IR photometry at 3.5, 8.4, and 11 jim. The 
color 3.5 jim - 8.4 jim was taken as a measure 
of silicate emission and the color 8.4 jim - 11 
jim as a measure of optical depth. Gehrz and 
Woolf avoided many of the problems of the 
data analysis by adopting a model of a thin 
spherical shell with temperature T at a distance 
R from the center of the star. When the color 
8.4 jun - 1 1 jim is plotted against the color 3.5 


fim - 8.4 jim, comparison with model curves 
gives both r( 11 jim) and ( R/R *) (77T*) _1/2 . 
Optical depths for class la supergiants were 
found to be 0. 1 to 0.2, as compared with about 
unity for Miras and about 0.25 for SRb vari- 
ables. The adoption of a dust temperature of 
900 K for all calculations led to estimates for 
the inner radii of shells of about 4 /?* for typical 
red giants and Miras and about 20 /?* for 
supergiants. 

In a study of mass loss from nine stars, 
which included giants, Miras, and supergiants, 
Hagen (1978) computed optically thin model 
dust shells (Dyck and Simon, 1975) to calculate 
infrared fluxes, making some allowance for 
self-absorption. After calculating a number of 
emission profiles as a function of r for assumed 
values of the other model parameters, Hagen 
chooses a value of r by comparing the predicted 
and observed emission bump heights and 
evaluates the dust column density from Equa- 
tion (5-8). The hydrogen column density then 
follows from the assumed grain composition. 
In calculating mass-loss rates from total column 
densities, Hagen (1978) took as inner shell radii 
the distances at which silicate grains condense, 
which vary from 1.2 R * for a 2500 K star to 
1 .9 R* for a 3600 K star. As Hagen points out, 
these small distances are inconsistent with the 
analysis of the optical line profiles, which were 
modeled on the assumption that the radius of 
the inner shell is larger than /?* by at least a 
factor of 3. In a more recent survey of 42 red 
giants and supergiants (Hagen et al., 1983), the 
inner shell radius is set at a uniform value of 
10 R* for all stars. 

In the same paper, Hagen (1978) also derives 
hydrogen column densities from optical lines, 
using Sr II for the cooler stars later than M5 
and Sr II, Ba II, and Sc II for the early M stars. 
Since she assumes that the same fraction of all 
elements is condensed into grains, the addition 
of the two values of 7V(H) gives the total col- 
umn density of CS matter, and the ratio of the 
two is the relative proportion of metals existing 
as gas and dust, respectively. It is hard to 
understand how two elements such as Sr and 
Si, with fractional abundances differing by five 
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orders of magnitude, can be condensed on 
grains in the same proportions. At the very 
least, the investigation should be repeated, us- 
ing optical observations of elements with abun- 
dances comparable to that of Si. 

One of the principal conclusions from these 
studies is that the total quantity of CS gas de- 
rived from the Sr II 4077 A line is uncorrelated 
with the quantity of CS dust implied by the 
strength of the 10-fim silicate feature. Such a 
correlation is expected from the theory of dust- 
driven mass loss (see Deguchi, 1980), but its 
absence may not be too significant in view of 
the many uncertainties in the derivation of dust 
column densities. Moreover, the use of a single 
line, 4077 A of Sr II, as a measure of total gas 
content implies far more confidence in ioniza- 
tion theory and in solar-like abundances than 
is warranted. It appears that the fraction of 
metals in the form of dust and gas, respective- 
ly, are not at all well determined. 

Submillimeter Dust Emission. In principle, the 
rate of mass loss in the form of dust may be 
estimated more directly from the thermal emis- 
sion of the dust at submillimeter wavelengths 
than from the 10-^m silicate feature. Sopka et 
al. (1984) have reported observations of 17 
evolved stars, obtained with the United King- 
dom (U.K.) infrared telescope on Mauna Kea 
at an effective wavelength of 400 /im. At this 
wavelength, the dust is optically thin, and 
moreover, the dust emission follows the 
Rayleigh-Jeans law. Consequently, uncertain- 
ties in grain temperature are not as serious as 
at shorter wavelengths. Nine of the objects were 
detected, and rates of mass loss, varying in 
amount from about 10~ 7 to 7.6 x 10' 6 
M 0 /yr, were estimated from the fluxes and are 
listed in Table 5-2. The greatest potential 
sources of error are the distance and *(400 ^m), 
the cross section per unit mass of the grains. 
Although the composition of both carbon and 
silicate grains is unknown, several lines of 
evidence (Draine, 1981; Rowan-Robinson and 
Harris, 1983b; Jura, 1983a) suggest that car- 
bon in the interstellar medium is amorphous 
rather than crystalline in structure and that the 


silicates are likely to be amorphous as well. The 
slopes of the far-IR to submillimeter spectra of 
both carbon-rich and oxygen-rich stars also 
support the conclusion that both carbon and 
silicate grains have amorphous structures. 
Thus, Sopka et al. (1984) adopt a common 
value of x — 20 cm 2 for both carbon-rich and 
oxygen-rich stars, which is estimated to be ac- 
curate to within a factor of 3, and is about an 
order of magnitude less for crystalline 
materials. The seventh column of Table 5-2 
gives the ratio of mass lost in the form of gas 
and dust, respectively, which is seen to be con- 
sistent with the interstellar value of 100. 

It is noteworthy that the largest values of 
M d are found for the bipolar objects, OH 
231.8 4 - 4.2, CRL 2688, CRL 618, the plane- 
tary nebula NGC 7027, and VY CMa, in all of 
which the flows are highly anisotropic. It might 
be believed that analysis of such objects based 
on spherical symmetry might give incorrect 
results, but Jura’s (1983b) calculation of mass 
loss from oblate ellipsoids shows that the rates 
are not greatly different from values derived by 
assuming spherical symmetry. 

Millimeter Waves 

Since it was first detected a few years ago 
(Zuckerman et al., 1977, 1978; Lo and Bechis, 
1977; Morris et al., 1979; Morris, 1980; see also 
Rieu, this volume, for an up-to-date review), 
thermal millimeter-wave radiation in rotational 
lines of molecules has become an increasingly 
important tool for studying mass loss from red 
giants. The CO lines have been most widely 
used in rate determinations. The 1-0 line has 
been detected in 50 stars, mostly of types M, 
S, and C (Knapp and Morris, 1984), and the 
2-1 line in 17 of the same objects (Knapp et al., 
1982). Such lines are particularly valuable for 
rate determinations because they are formed far 
from the star, where the flow is likely to have 
attained both its terminal velocity and spherical 
symmetry. One disadvantage is that they appear 
in emission, and hence the mass of CO in the 
envelope must be derived from detailed models 
in which the square of the distance is a free 
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Table 5-2 

Mass-Loss Rates and Expansion Velocities of Selected Red-Giant Stars 


Star 

Spectral 

Class 

M 

9 

\/ 

exp 

(km/s) 


M //W w 

9 d 

Ref. t 

q Per 

M4 11-111 

0.012 

7 



2 

1 1 9 Tau 

M2 lb 

0.24 

7 



2 

a Ori 

M2 lab 

4 

10 



3,4, 5 

6 Gem 

M 1 -M2 la-lab 

2.1 

8 



2 

a Sco 

Ml. 5 lab 

5.4 

9 



2,6 

a Her 

M5 lb— II 

0.11 

8 



2,7 

6 2 Lyr 

M4 II 

0.048 

7 



2 

R Lyr 

M4-5 III 

0.014 

5 



2 

fi Cep 

M2 la 

1.0 

9 



2,4 

6 Sge 

M2 II 

0.02 

28 



8 

IRC + 10011 

M8 

14 

24 

0.07 

200 

9,10 

o Cet 

M6e 

0.65 

4 



9 

NML Tau 

M 

5.1 

28 



9 

S CMi 

M7e 

6.3 

18 



1 

RS Cnc 

M6 lb— II 

0.29 

11 



9 

R LMi 

M7e 

1.0 

6 



1 

R Leo 

M8e 

0.092 

7 



9 

RX Boo 

M7e 

0.38 

8 



9 

NML Cyg 

M6 III 

1.8 

21 

0.058 

31 

1,10 

R Cas 

M7e 

0.66 

11 



1 

R And 

S6e 

0.37 

15 



9 

W Aql 

S 

9.8 

20 



9 

X Cyg 

S7e 

0.18 

8 



1 

1C 418 

PN 

5.0 

17 



1 

NGC 6543 

PN 

5.1 

14 



1 

NGC 7027 

PN 

110 

22 

1.3 

85 

9,10 

R Scl 

C 

31 

16 



9 

IRC + 50096 

C 

6.3 

18 



9 

CRL 618 

C 

77 

18 

0.85 

90 

9,10 

CRL 865 

c 

23 

14 



9 

IRC + 10216 

C9 

55 

17 

0.17 

320 

9,10 

CIT 6 

C4 

3.0 

17 



9 

IRC - 10236 

C 

47 

10 



9 

V Hya 

C6 3 e 

4.6 

18 



9 

IRC + 20326 

C 

23 

10 



9 

CRL 2155 

c 

17 

20 



9 

CRL 2199 

c 

13 

15 



9 

IRC + 20370 

C7 

10 

9 



9 

V Cyg 

C7e 

25 

14 



1 

CRL 2688 

C 

160 

19 

2.3 

70 

9,10 

IRC + 40485 

C 

17 

13 



1 

CRL 3068 

C 

7.3 

14 

0.086 

85 

9 

IRC + 40540 

C8 

24 

15 



9 

OH231.8 + 4.2 


130 

68 

7.6 

17 

9,10 


* m and M are rates derived for gas and dust, respectively, in units of 10 M 0 I yr. 

T References' 1 . Knapp et al. (1982); 2. Sanner (1976); 3. Bernat et al. (1979); 4. Hagen et al. (1983); 5. Mauron et 
al. (1984); 6. Hjellming and Newell (1983); 7. Reimers (1977); 8. Reimers and Schroder (1983); 9. Knapp and Mor- 
ris (1984); 10. Sopka et al. (1984). 
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parameter; another is the uncertainty in the 
composition ratio CO/H 2 , which may, in 
some cases, be derived indirectly from observa- 
tion (see below). Circumstellar chemistry is a 
new field of research, which is not yet on a firm 
quantitative basis. (See the review by Glassgold 
and Huggins in this volume.) 

Figure 5-6 shows typical spectra in the region 
of the CO J = 2-1 line for Betelgeuse and 
IRC + 1001 1 . As shown by Morris (1975), the 
profiles of emission lines from unresolved 
spherical gas clouds expanding at constant ve- 
locity are flat-topped, as in Betelgeuse, when 
the cloud is optically thin, and parabolic, as in 
IRC + 10011, when the cloud is optically thick. 
The width of the profile at zero intensity or 
power is twice the expansion velocity, K , 
whereas the velocity at the center of the line is 
the systemic, or center-of-mass velocity of the 
star, k. For Betelgeuse, the derived velocities 
are V ex P ~ *5.3 ± 2.5 km/s and K = 18.8 ± 
2.5 km/s (Knapp et al., 1980), in good agree- 
ment with the optical values V = 16.5 and 

ex P 

10.0 km/s (Goldberg et al., 1975; Bernat et al., 
1979) and K = 21 km/s (Adams, 1956). The 
expansion velocity derived from the emission 
profile apparently corresponds to the larger of 
the two displacements observed in optical ab- 
sorption spectra of a Ori. The absence of emis- 
sion from the inner shell can be understood 
(Knapp et al., 1980) if the shell is much closer 
to the star than proposed by Bernat et al. 
(1979), 8 rather than 150 R *. Moreover, 
spatial mapping of the millimeter-wave emis- 
sion from Betelgeuse gives a radius of about 400 
R* for the outer shell, again much smaller than 
the value of about 2000 R * found by Bernat 
et al. (1979). Goldreich (1980) has pointed out 
that the gas in the circumstellar shell is heated 
by the drag force between gas and dust and not 
by the thermal dust temperature. It is therefore 
incorrect to equate the excitation temperature 
with the equilibrium temperature of the dust, 
at least for those objects in which the rate of 
mass loss and therefore the density are relatively 
low (Morris, 1980). 

Theoretical models of circumstellar enve- 
lopes around red giants, in which rates of mass 



Figure 5-6. Spectra in the region of the CO 
J = 2-1 line at 230 GHz for a Ori (top) and 
IRC + 10011. Possible blends with lines of 
other molecules are indicated (Knapp et al., 
1982). 


loss and chemical compositions are parameters 
to be derived from molecular line observations, 
have been calculated by Goldreich and Scoville 
(1976) for OH/IR stars, by Kwan and Hill 
(1977) and Kwan and Linke (1982) for the car- 
bon star, IRC + 10216, and by Morris (1980; 
see also Morris and Alcock, 1977) for a wide 
range of stellar parameters. In all of these 
models, mass loss is assumed to be driven by 
radiation pressure on dust grains in a spher- 
ically symmetric atmosphere. The gas temper- 
ature distribution is determined by the balance 
between heating by collisions with grains and 
cooling by adiabatic expansion and by the 
emission of H 2 0 in O-rich stars and of CO in 
carbon stars. Because most of the CO emission 
occurs hundreds of stellar radii from the star, 
Kwan and Linke (1982) assume constant out- 
flow velocity and use the Sobolev approxima- 
tion to solve the radiative-transfer problem. 
The solution contains three free parameters: M, 
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the rate of mass loss, Q, the momentum 
transfer efficiency factor averaged over the con- 
tinuum energy distribution, and (CO)/(H 2 ), 
the abundance ratio carbon monoxide to mo- 
lecular hydrogen. It is expected that, in the cir- 
cumstellar envelopes of the cooler late-type 
stars, hydrogen will be predominantly in mo- 
lecular form (Knapp and Bowers, 1983). The 
dust grains are assumed to constitute 1 percent 
of the mass and to have a density of 1 gm/cc. 
Since these quantities always enter with Q, their 
uncertainties are lumped together with the un- 
certainty in Q , which may be checked by cal- 
culating the optical depth of the dust and com- 
paring it with that of infrared observations. 

Morris (1980; see also Morris and Alcock, 
1977) bases his model on temperature and ve- 
locity distributions derived from their dust 
model by Goldreich and Scoville (1976). For ex- 
ample, the velocity law is: 

V * F - 0.37 X 10 16 /r , (5-9) 

lerm 

for r > 10 15 cm, equivalent to 20 /?* for 
Betelgeuse. 

Morris solves the coupled equations of radi- 
ative transfer and statistical equilibrium to cal- 
culate the radial distribution of the populations 
in the lowest 1 1 rotational levels of both the 
first and second vibrational levels of CO. The 
processes contributing to the level excitations 
are: 

1. The millimeter-wave radiation field, in 
part from the 3 K background and in 
part from diffused CO line radiation. 

2. The 4.6-/im radiation field exciting the 

v = 0-1 transitions. The central star is 

assumed to be a blackbody with a tem- 

perature of 2000 K and a radius of 5 x 
10 13 cm, and a parameter, W, is sup- 

posed to measure the amplification of 

the 4.6-fim radiation by the surrounding 

dust. No account appears to have been 

taken of the weakening of the stellar ra- 
diation field by CO absorption in the 
photosphere. 


3. The levels are also excited by collisions 
with H 2 molecules, the results being in- 
sensitive to the assumed temperature 
distribution. 

Results are presented for the central inten- 
sities of the transitions J = 1-0, 2-1, and 3-2 
as a function of M (the mass-loss rate), W (the 
measure of 4.6-^m luminosity), and/(the frac- 
tional abundance of CO relative to H 2 ), and 
for the beamwidths of the 1 1-m National Radio 
Astronomy Observatory telescope and the 10-m 
Owens Valley Radio Observatory telescope. For 
a completely unresolved envelope, the bright- 
ness temperature is inversely proportional to the 
distance squared. Actually, both of the above 
telescopes can resolve envelopes at distances 
closer than 1000 parsecs, which means that an 
increase in distance is partially compensated by 
an increase in the surface area intercepted by 
the beam. Therefore, the peak surface bright- 
ness falls off less rapidly than as D 1 . 

The results show that, when the product/ 
x M < 10 -9 , the envelope is optically thin 
both in the rotational lines and in the 4.6-^im 
lines. The line profiles are flat for D >1 kpc 
and become double-peaked as D decreases be- 
cause the outer parts of the emitting region are 
no longer within the telescope beam. When/ 
x M > 2 x 10" 2 * * * * * * 9 , the envelope is optically 
thick in both rotational and vibrational lines. 
Increased density causes the excitation to be 
dominated by collisions with H 2 . Finally, 
when / x M ~ 10~ 9 , weak maser emission 
may occur. 

It might be believed that only the product 
/ x M could be derived from comparing model 
calculations with observations, but that is 
usually true only in cases in which the levels are 
excited purely by radiation. Otherwise, the 
same value of the product / x Ki may produce 
profiles of quite different intensity and shape 
(double-peaked, flat, or parabolic) depending 
on the relative values of / and Ki. The reason 
is that the value of M is a measure of the den- 
sity, which controls the rate of collisions and 
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thus the excitation temperature. An increase in 
excitation temperature increases the size of the 
observed emitting region and therefore modifies 
the shape of the profile. Thus, detailed model- 
ing of emission-line profiles can yield both the 
mass-loss rate and information on the chemical 
composition. 

Another type of modeling, which seeks to 
combine optical and radio observations, was 
conceived by Jura and Morris (1981) and ap- 
plied to the shell of Betelgeuse. The model is 
derived from observations of the 2-1 line of CO 
(Knapp et al., 1980) combined with measures 
of the surface brightness of the shell in the 
monochromatic radiation 7699 A of K I at 
distances of 20 to 60 arc-sec from the star 
(Honeycutt et al., 1980). At such great distances 
from the star, the ionization of K is controlled 
entirely by the interstellar radiation field, al- 
though modified by dust extinction, or shield- 
ing, in the shell itself, and therefore the uncer- 
tainty about the stellar radiation field is 
avoided. The analysis has recently been re- 
peated by Mauron et al. (1984) using new sur- 
face brightness measurements made with a 
Fabry-Perot etalon and a charge-coupled device 
(CCD) camera. The method yields the radial 
density distribution of potassium and of hydro- 
gen if the solar abundance ratio is assumed. A 
by-product of the analysis is the abundance 
ratio carbon/potassium. So far, such optical 
surface brightness distributions have been 
measured only for the star Betelgeuse. 

In their survey of mass loss from 17 evolved 
stars, Knapp et al. (1982) avoided detailed 
modeling by putting the stars into two cate- 
gories, those with optically thick and optically 
thin envelopes, respectively, and deriving sim- 
ple expressions for the mass-loss rate in terms 
of antenna temperature, T * A , expansion veloc- 
ity, V, and distance, D. The expressions are 
designed to fit the IRC + 10216 model of Kwan 
and Hill (1977) and Kwan and Linke (1982) as 
a prototype for the optically thick envelopes 
and the a Ori models of Bernat et al. (1979) and 
Jura and Morris (1981) as typical for optically 
thin envelopes. The models predict that the 
ratio of antenna temperatures of the 2-1 and 


1-0 lines of CO should be about 3 in optically 
thick envelopes and several times larger for op- 
tically thin shells. According to this criterion, 
15 of the 17 shells were shown to be optically 
thick, the only exceptions being Betelgeuse and 
Mira. 

In their optically thick approximation, 
Knapp et al. (1982) assume that each source is 
a sphere of constant surface brightness, the 
radius being fixed by the gas density at which 
the CO lines are no longer thermalized. How- 
ever, Jura (1983b) has shown that the outer op- 
tically thin regions of the sources cannot be 
neglected, since they contribute substantially to 
the total flux. Allowance for this effect reduces 
Knapp et al.’s mass-loss rates by a factor of 2 
to 3. A further reduction of about 30 percent 
becomes necessary when the effects of dust 
shielding and molecular self-shielding are con- 
sidered (Jura, 1983b). Jura’s analysis leads to 
a less steep dependence of mass loss on dis- 
tance, Z) 3/2 , rather than D 2 as found by Knapp 
et al. (1982). 

A new survey by Knapp and Morris (1984) 
has detected 50 sources in the 1-0 line. The 
analysis has been carried out with a minimum 
of simplifying approximations, and the results 
are probably more reliable than those of Knapp 
et al. (1982). As before, general expressions for 
T* a were derived for both optically thin and 
optically thick cases, respectively. For both 
cases, the antenna temperature is a function of 
M , Z), /, and r , the maximum extent of CO 
in the envelope. In optically thin objects, T* A 
also depends on the IR flux, W. For most 
objects, the model profiles are in satisfactory 
agreement with observations, but a few opti- 
cally thin shells, including those of o Ceti and 
a Ori, presented problems that require further 
study. In view of the many free parameters in 
the theory, the results must still be regarded 
with some caution. 

Microwave Emission 

VLA measurements of radio emission at 
centimeter wavelengths from the supergiants 
Betelgeuse (Newell and Hjellming, 1982) and 
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Antares (Hjellming and Newell, 1983) are con- 
sistent with emission of an optically thick ther- 
mal radio spectrum emitted by an extended 
slightly ionized chromosphere extending to 
several stellar radii. In addition, the observa- 
tions suggest that the wind from Antares is 
ionized by the Lyman continuum radiation of 
the companion B2.5 V star. Figures 5-7a and 
5-7b are 4.885-GHz images of the Sco system 
made with the VLA, showing the unresolved 
M star on the left and the nebulosity sur- 
rounding the B star, the location of which is 
marked by a cross. Figure 5-7c is a theoretical 
model of the ionized portion of the Antares 
stellar wind, which is created by the ultraviolet 
flux from the B star. The agreement between 
theory and observation is most convincing and 
suggests a new method for measuring the rate 
of mass loss from an M supergiant. Hjellming 
and Newell (1983) point out that, at the ioniza- 
tion front between the neutral and ionized hy- 
drogen, the number of neutral hydrogen atoms 
flowing per square centimeter across the ioniza- 
tion front must equal the attenuated Lyman 
continuum photon flux which ionizes the gas. 
Thus, 

N l v i = L(ly) • exp(-T,)/4ir • r, 2 , (5-10) 

where N j and v ^ are the hydrogen density and 
velocity of the wind at distances of R ] and r { 
from the M and B stars, respectively; t, is the 
optical depth for the Lyman radiation between 
the ionization front and the star, assumed to 
be unity; and L(ly) is derived from a model of 
the microwave radio emission from the star. Its 
value is proportional to the square of the stellar 
distance and to the -1/6 power of the electron 
temperature. Since the mass-loss rate is: 

M = 4ir ■ e, ■ ”, ' > (5-1 1) 


where q is the density, the expansion velocity 
may be eliminated, and M is given by: 


to = (e/Wj) • L(ly) • exp(-T,) 

• ( V r ,) 2 • 

Thus, the mass-loss rate may be derived entirely 
from observed quantities and straightforward 
application of the theory of H II regions. 

Because systems like a Sco are uncommon, 
Hjellming and Newell’s (1983) method will not 
have wide application. Moreover, only the hot- 
test late-type giants will have winds in which 
most of the hydrogen is ionized. Thus, the 
microwave flux from M giants and supergiants 
will not be useful for the determination of 
mass-loss rates unless there is an independent 
means for measuring the fraction of ionized 
hydrogen, a Ori is known to emit weak 
microwave radiation (Newell and Hjellming, 
1982) from an extended warm chromosphere 
in which the flow is probably accelerating, but 
in which the fractional ionization cannot be 
more than 0.01 to 0.10. The same fraction has 
been found for a 1 Her (Drake and Linsky, 
1984). Spergel et al. (1983) observed 31 evolved 
red giants with known high rates of mass loss 
with the VLA at 5 GHz and detected emission 
from only four stars. Upper limits of about 
10~ 3 were derived for the fractional hydrogen 
ionization. On the other hand, the probable 
detection of a Boo, also with the VLA at 6 cm 
(Drake and Linsky, 1984), implies a mass-loss 
rate of 1 x 10“ 10 Af © /yr, which indicates that 
hydrogen in the wind is well ionized. 

OH/IR Masers 

The subject of red-giant masers has recent- 
ly been reviewed by Elitzur (1981). Maser radia- 
tion at radio frequencies, which is usually found 
in oxygen-rich giants and supergiants of spec- 
tral types later than M5, is emitted by SiO, 
H 2 0, and OH molecules. The SiO radiation is 
emitted by excited vibrational levels, which are 
pumped by collisions, and therefore must be 
located in the lower atmosphere of a star. The 
OH lines, at 1612, 1665, 1667, and 1720 MHz, 
originate from the ground vibrational state and 
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Figure 5-7. Images of the a Sco radio source made with the VLA at 4.885 GHz. The unresolved M 
star is at the left, and the position of the B star is marked by a cross at the right, (a) and (b) are 
images constructed with different methods of weighting, and (c) is a theoretical model of an H II 
region produced in the wind of the M star by UV photons from the B star (Hjellming and Newell, 1983). 


are pumped by IR radiation, as evidenced by 
the variations in phase of IR and OH emission 
(Harvey et al., 1974). Three of these four lines 
have been found to emit maser radiation in late- 
type stars, the so-called main lines at 1665 and 
1667 MHz, and the satellite line at 1612 MHz. 
Those stars with stronger main lines are clas- 
sified as Type I and those with stronger satellite 
lines as Type II. We shall be concerned here 
only with Type II masers because model cal- 
culations show that the satellite line is emitted 
only by stars with very high rates of mass loss. 
The enormous extent of maser shells has been 


verified independently by different methods. 
Thus, Baud (1981) observed, with the VLA, ex- 
tended OH emission from the source OH 26.5 
+ 0.6 over a diameter of 4.6 x 10 16 cm. 
Jewell et al. (1980), making use of the phase 
relation between variations in the OH 1612- 
MHz emission originating from the front and 
back of the shell, derived a diameter of (6.6 ± 
1.4) x 10 16 cm for the maser shell around IRC 
+ 10011 . 

Some typical spectra in the region of the 
1612-MHz line are shown in Figure 5-8. The 
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two peaks originate from the approaching and 
receding halves of the shell, respectively. 
Hence, the separation of the peaks measures 
twice the expansion velocity, and the velocity 
half way between the peaks is the systemic 
velocity of the star. Mass-loss rates are deter- 
mined from these observations on the assump- 
tion that the mechanism is radiation pressure 
on dust grains, which is highly likely for stars 
with such thick dust shells. If radiation pressure 
on grains is much larger than all other forces, 
such as gravitation, it can be shown that the 
equation of motion, when combined with the 
equation of continuity and the definition of op- 
tical depth (Elitzur, 1981), reduces to: 


M = L • t d /Vc , (5-13) 


where L is the mean IR flux, V is the terminal 
expansion velocity, c is the velocity of light, and 
t d is the effective optical depth of the dust 
through the shell in response to radiation 
pressure. Its value is on the order of 2 for dense 
shells. 

A second expression for the mass-loss rate 
may be obtained by integrating the equation of 
motion to obtain an expression for the terminal 
velocity, which, when inserted into the equa- 
tion of continuity, gives (Salpeter, 1974): 



VELOCITY (km/s) 


Figure 5-8. Spectra of Type II OH sources at 
1612 MHz: Top: OH 1.5-0.0; middle: OH 
2.2-1 . 7; bottom: OH 2. 6-0. 5. The left and right 
peaks arise from the approaching and receding 
halves of the shell , respectively (Baud et ai, 
1979). 


M = 2tt • R o - V • t d /x d • / , (5-14) 

where t d is the optical depth at a given 
wavelength, say, 9.7 /an, x is the mass absorp- 
tion coefficient at the same wavelength, /is the 
fraction of the total mass in the form of dust, 
and R is the radius at which dust forms. For- 

o 

rest et al. (1978) have used both equations to 
obtain estimates of mass loss from OH26.5 
+ 0.6. From Equation (5-14), with r>6, *(9.7 
ftm) = 3000, RJr* >3, and r* = 6.5 x 10 13 
cm, M>1.5 x 10~ 5 M©/yr. Similarly, Equa- 
tion (5-13), with L = 10 4 , which is appropriate 
for Miras, V = 13 km/s, and r = 3, gives M 


>4.5 x 10‘ 5 A/©/yr. Mass-loss rates 
calculated for OH/IR stars with Equation 
(5-13) appear in Table 5-3. 

Results 

Compilations of mass-loss rates for late-type 
stars, derived by the methods described above, 
will be found in Gehrz and Woolf (1971), 
Reimers (1975a, 1975b), Sanner (1976), Bernat 
(1977), Hagen (1978), Werner et al. (1980), 
Knapp et al. (1982), Hagen et al. (1983), Jura 
(1983b), Baud and Habing (1983), Engels et al. 
(1983), Sopka et al. (1984), Pottasch (1984), 
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Table 5-3 

Periods, Expansion Velocities, and Mass-Loss Rates of OH/IR Stars 
Calculated from Equation {5-13) 


Star 

Period 

(days) 

\/ 

exp 

(km s ) 

M 

HO' 6 /W 0 /yr) 

Ref.* 

R Aql 

280 

6.5 

0.8 

1 

U Ori 

372 

4 

0.3 

1 

RR Aql 

394 

6.5 

0.5 

1 

WX Ser 

425 

7.5 

1 

1 

OH31.7 - 0.8 

510 

13 

10.9 

2 

OH26.4 - 1.9 

540 

12 

5.2 

2 

OH 1 9.2 - 1.0 

610 

15 

24 

2 

OH28.7 - 0.6 

640 

18 

22.4 

2 

OH45.5 + 0.1 

720 

18 

4 

2,3 

OH39.9 + 0.0 

770 

15 

28.7 

2 

OH30.1 -0.2 

970 

18 

23 

2,3 

0H30.7 + 0.4 

1 140 

18 

8.1 

2 

OH26.2 - 0.6 

1330 

22 

33.2 

2 

OH26.5 + 0.6 

1630 

14 

75 

2,3 

OH32.8 - 0.3 

1750 

16 

60 

2,3 


•References: 1. Pottasch (1984); 2. Engels et al. (1983); 3. Werner et al. (1980). 


and Knapp and Morris (1984), among others. 
A useful tabulation of circumstellar line 
strengths and velocity displacements for 61 red 
giants is found in Boesgaard and Hagen (1979). 
Even for such well-studied stars as a Ori, a Sco, 
and IRC + 10216, published values of the 
mass-loss rates differ by factors of 10 to 100, 
and hence, it does not seem useful to present 
them all in tabular form. In many cases, how- 
ever, discordant values can be weeded out for 
known reasons, and the remaining values, al- 
though still uncertain by factors of perhaps 2 
to 10, serve a useful purpose in setting con- 
straints on mechanisms and evolutionary 
processes. 

Thus, Tables 5-2 and 5-3 are compilations 
of mass-loss rates of evolved stars selected, 


because they seemed to be relatively well deter- 
mined, from among several different classes of 
stars. Table 5-2 contains expansion velocities 
and mass-loss rates for early- and late-type M 
stars, for carbon stars, and for a few S stars 
and planetary nebulae. Rates for the early M 
stars have been derived mainly from optical 
data, and the selection of the best values can 
be a matter of personal taste, a Ori is a case 
in point. Jura and Morris (1981) and Mauron 
et al. (1984) used exactly the same method of 
analysis, combining data on CO line emission 
with surface brightness measurements in neutral 
potassium, but their rates differ by a factor of 
4 because the surface brightness distribution of 
K I emission measured by Mauron et al. (1984) 
differed considerably from that of Honeycutt 
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et al. (1980). It should be noted that the abun- 
dance ratio K/C found by Jura and Morris 
(1981) was 25 times larger than the solar value, 
but Mauron et al. (1984) reduced this discrep- 
ancy by a factor of 6. A reduction in the C/H 
ratio by a factor of 1 to 5 would not be unex- 
pected for processed material. For the later 
spectral types and the carbon stars, the most 
reliable rates are probably those from rotational 
lines of CO. When available, the most recent 
values of Knapp and Morris (1984) have been 
adopted. The data show several clear trends: 

1. Among the early M stars, rates for the 
supergiants exceed those for the giants by 
a factor of 10 or more. 

2. Rates for the Miras and planetary 
nebulae are comparable to those for the 
early-type supergiants, but the value for 
NGC 7027 is at the extreme upper end 
of the range. 

3. Knapp et al. (1982) found that the rates 
for the carbon stars were 25 times larger, 
on the average, than those for the O-rich 
stars, but this result is not borne out by 
the newer analysis of Knapp and Morris 
(1984), who allowed for the presumably 
higher CO abundance in C-rich stars and 
included the effects of excitation by stel- 
lar IR on the optically thin objects. Thus, 
if/ = 3 x 10"* for oxygen Miras and 8 
x 10" 4 for the carbon stars, the range of 
rates of mass loss found for the two 
groups of stars is about the same. The 
results of the Knapp and Morris (1984) 
survey are summarized in their Figure 16 
in the form of histograms in which the 
number of stars per pc 3 is plotted 
against Kf. 

4. The M and S stars show considerable 
scatter in their expansion velocities, with 
values from 4 to 28 km/s, whereas 
among the carbon stars, u exp is essen- 
tially confined between 10 and 20 km/s. 


Table 5-3 is a compilation of periods, expan- 
sion velocities, and mass-loss rates of OH/IR 
stars, the rates having been calculated from 
Equation (5-13). The table contains four 
OH/IR objects optically identified as Miras, 
together with 1 1 unidentified OH/IR stars. The 
results for the OH/IR stars can be better under- 
stood in view of recent studies of these fasci- 
nating objects. About 350 OH masers are 
known (Engels, et al., 1981), of which less than 
one-third have been identified optically, most 
as M-type Miras but also M supergiants, semi- 
regular variables, and some very red IRC 
sources. A recent survey of about 200 OH/IR 
stars (Engels et al., 1981) revealed that all 
OH/IR stars are variable, with amplitudes up 
to two magnitudes in the K band. At least one- 
half the stars have longer periods than those of 
Miras or IRC sources, with values up to 1600 
days or more. It now seems to be well estab- 
lished (Baud and Habing, 1983; Engels et al., 
1983) that the optically identified OH/IR stars 
and the unidentified sources form a continuous 
sequence of Mira variables extending from the 
optical objects with periods up to about 400 
days to the unidentified sources with periods 
up to nearly 2000 days. As seen in Table 5-3, 
the expansion velocities along this sequence in- 
crease on the average from <10 km/s for opti- 
cal Miras to 22 km/s. From statistical studies 
of velocity dispersions and galactic concentra- 
tions, Baud et al. (1981) found the expansion 
velocities to be correlated with age and main- 
sequence mass in the sense that small velocities 
go with small mass. 

Thus, the progression with period is also one 
with main-sequence mass from about 0.5 M© 
for the optical Miras to 6.5 M© for the OH/IR 
sources with the longest periods. The mass-loss 
rates from unidentified OH/IR stars are also 
strikingly higher, and the increase occurs rather 
abruptly at periods over 400 days. There ap- 
pears to be a threshold period beyond which 
the rate of mass loss suddenly increases by a 
factor of 10, and the star becomes enveloped 
in an opaque dust cloud. 
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NONSTEADY MASS LOSS 
Selected Examples 

Evidence is accumulating that, in some stars 
at least, mass loss may occur sporadically or 
at a variable rate. The following examples are 
illustrative rather than complete. 

Multiple components of CS lines have been 
observed in the spectra of a number of red gi- 
ants and supergiants, chiefly the latter (Sanner, 
1976; Bernat, 1981). From a study of CO lines 
in the 4.6-jim band in nine red giants and 
supergiants, Bernat (1981) concluded that the 
multiple components originated in separate 
mass-loss episodes rather than as intermittent 
modulations of continuous mass loss, which is 
favored by Ridgway (1981). The matter is still 
controversial. 

Reimers (1975a, 1975b) noted variability in 
the K-line profiles of early M giants, stars lying 
near the transition zone in the HR diagram be- 
tween stars with low and high rates of mass 
loss. 

Brooke et al. (1974) found evidence for 
large-scale motions in the atmosphere of a Ori, 
which appeared to be sporadic in time or posi- 
tion on the stellar surface. 

HR 8752 and q Cas, two F la supergiants 
with extensive circumstellar shells containing 
little dust, have undergone dramatic changes in 
their atomic and molecular spectra in the 
2.3-fim region over a period of years (Lambert 
et al., 1981). During a 5-year period from 1975 
to 1979, P Cygni profiles in HR 8752 were 
transformed into inverse P Cygni profiles, sug- 
gesting the fallback of a shell ejected in 1975. 
The mass-loss rate derived for q Cas is very 
large, exceeding 1CT 3 M 0 /yr, leading Lambert 
et al. (1981) to propose that it is an outburst 
connected to a peculiar decrease in brightness 
observed in 1946/1947. 

If the stationary atmospheric layer observed 
by Hinkle et al. (1982) in x Cyg and other Miras 
is the source from which mass escapes, the vari- 


ations in the thickness of this layer in \ Cygni 
may be an indication of variable mass loss. 

The far-infrared (40 to 250 fim) size of IRC 
+ 10216 (Fazio et al., 1980) suggests more dust 
at large distances than would be expected for 
an r' 2 density distribution. This implies that, 
about 2000 years ago, the rate of mass loss was 
twice its present value. Mitchell and Robinson 
(1980) also derived a relatively shallow density 
gradient, varying as r -13 . On the other hand, 
the brightness distribution at 1 1 /*m is consis- 
tent with r~ 2 (Sutton et al., 1979). Rowan- 
Robinson and Harris (1983b) conclude, 
however, that errors of measurement are too 
large to rule out the r' 2 density law. 

Optical spectropolarimetry of dust-en- 
shrouded extreme carbon stars reveals that their 
circumstellar shells are highly organized (Cohen 
and Schmidt, 1982). In particular, GL 1403 = 
IRC + 30219 = CIT 6 seems to show a torus 
of cool dust and hot bipolar lobes, replenished 
by outbursts occurring after the minima of the 
640-day light variation. Particularly striking 
were the spectral changes in 1977 through 1980 
(Cohen, 1980), when the stellar photosphere 
nearly disappeared and the spectrum consisted 
of emission lines of H, He, and forbidden lines. 
Cohen (1980) suggests that these lines were 
formed when a recent episode of mass loss col- 
lided with the existing circumstellar envelope. 
These and similar objects may represent the 
evolution of high-mass carbon variables into 
bipolar nebulae. 

VLA observations of extended OH maser 
emission for the envelope around OH26.5 + 
0.6 show large-scale spatial asymmetries and 
perhaps a ring-like structure seen edge-on, in- 
dicating real density variations in the envelope. 

Approximately twelve Type II OH/IR 
sources have been found to have dust shells of 
high optical depth and to display deep silicate 
absorption at 10 ^m. Werner et al. (1980) have 
modeled these shells with an r~ l 5 density law, 
whereas Rowan-Robinson (1982) has proposed 
a double-shell model in which a cool star with 
a spherically symmetric dust shell expanding at 
constant velocity is surrounded by a second 
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shell of cold absorbing dust. The second shell 
is believed to represent an earlier episode of 
ejection. The mass of the cold dust is estimated 
as 4 to 15 Mq , which implies that the sources 
are very massive highly evolved stars. 

Near-infrared (1.65 to 4.8 fim) speckle in- 
terferometry of IRC + 10216 (Dyck et al., 
1984) shows the outer regions of the object to 
be circularly symmetric on scales > 1.6 inch, 
but at intermediate distances, the brightness dis- 
tribution is strongly asymmetrical, being more 
extended north-south than east- west, as found 
previously by McCarthy et al. (1980). McCar- 
thy et al. also reported variations with time of 
the apparent envelope size, which have been 
confirmed by Dyck et al. (1984). 

Evidence for Sporadic Mass Loss 
in ot Orionis 

a Ori has been relatively well observed for 
many years, and there is now an impressive 
amount of evidence for nonsteady ejections of 
matter at all levels of its highly extended at- 
mosphere and circumstellar shell. At the level 
of the photosphere, both visible light and radial 
velocity vary on the average in a period of 5.78 
years, apparently arising from pulsation, but 
these smooth variations are frequently dom- 
inated by irregular fluctuations, on a time scale 
of about 100 days, which often mask the varia- 
tions associated with the pulsation. (For a sum- 
mary, see Goldberg (1979, 1984).) Ordered 
changes in the linear polarization of the visible 
light have been observed over a number of years 
(Hayes, 1980, 1981, 1984; Tinbergen et al., 
1981). After a lapse of nearly 50 years, 
systematic monitoring of the photoelectric 
brightness of the star was resumed in 1979 
(Krisciunas, 1982; Guinan, 1983). As shown in 
Figure 5-9, reproduced from Guinan (1983), the 
amount of linear polarization and the blue 
magnitude are at least roughly correlated, 
which is consistent with asymmetric ejection of 
matter. Similarly, the absence of a one-to-one 
correlation between the short-term brightness 
and velocity variations (Sanford, 1933) makes 


it apparent that these changes are not global in 
nature. The wavelength dependence and time 
variation of the polarization are qualitatively 
consistent with the growth of a hotspot, or large 
convection cell, occupying a few percent of the 
area of the stellar disk (Schwarz and Clarke, 
1984; Clarke and Schwarz, 1984). 

The brightness distribution of emission from 
the inner shell of Betelgeuse within 10 stellar 
radii also shows pronounced asymmetries 
(Ricort et al., 1981; Goldberg et al., 1981; Rod- 
dier and Roddier, 1983) both in Ha light and 
in the continuum. The asymmetries persist to 
very large distance scales. Diffraction-limited 
images of the a Ori dust shell at 10 /xm, ob- 
tained with the NASA IR telescope on Hawaii 
(Bloemhof et al., 1984) show a marked asym- 
metry in the dust distribution, with peak inten- 
sity at a distance of 0.9 arc-sec from the star. 

The possible association of variability in 
brightness and radial velocity with mass loss 
suggests a look at the past behavior of a Ori 
for evidence of events that may have signaled 
the ejection of substantial amounts of matter 
from the photosphere (Goldberg, 1984). Obser- 
vations of visual brightness and radial velocity 
for this star have been carried out with vary- 
ing degrees of regularity and accuracy for well 
over 60 years, the most accurate and homo- 
geneous set of data having been acquired in 
1917 through 1931 (Stebbins, 1931; Sanford, 
1933). Although short-term changes in radial 
velocity and brightness were found to be un- 
correlated in detail, Sanford (1933) noted one 
or two major disturbances in which decreasing 
velocity was accompanied by increasing light 
and vice versa. 

The best example of such an event took 
place near the time of minimum radial veloci- 
ty in the 6-year cycle, at about Julian day 
2424500 (December 15, 1925). As shown in Fig- 
ure 5-10, a large decrease in radial velocity after 
day 4500 was accompanied by an initial increase 
in brightness, followed by rapid fading. It is as 
though gaseous matter were ejected from the 
star, then diffused outward and condensed in- 
to grains which became optically thick at visible 
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Figure 5-9 . B-magnitudes and percent polarization of light from Betelgeuse. Large open circles are 
seasonal means , and solid curve is Stebbins* (1931) mean blue light curve ( Guinan , 1983). 
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Figure 5-10 . Light and velocity variations of a 
Ori around the time of the unusual event of 
December 15, 1925, when a period of decreas- 
ing radial velocity was followed by a steep 
decline in visual brightness ( Goldberg , 1984). 


wavelengths. Similar dust-forming episodes ap- 
parently occur during the development of slow 
novae of the DQ Herculis type (Gehrz et al., 
1980) in which the visual extinction is accompa- 
nied by strongly enhanced infrared emission. 
Although photometric measures were not avail- 
able, the records of the American Association 
of Variable Star Observers strongly suggest that 
in November 1944 (Adams, 1956) and again in 
1978/1979 (Goldberg, 1979), similar events may 
have occurred at or soon after the minimum of 
the mean velocity curve. 

At radial velocity minimum, a pulsating star 
has reached its maximum velocity of expansion 
and is beginning to decelerate. Given the enor- 
mous extension of a supergiant atmosphere, its 
motion probably lags that of the star, and the 
atmosphere will tend to continue accelerating 
outward past the minimum, as has been ob- 
served. At this time, instabilities in the atmo- 
sphere may trigger mass ejection, followed by 
the condensation of dust grains. In the future, 
systematic monitoring of a Ori with a variety 
of techniques before, during, and after the 
velocity minimum may clarify the relationship 
between the short-term brightness and veloc- 
ity changes and the mass-loss process. 
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MECHANISMS 

Any acceptable theory of mass-loss mecha- 
nisms must explain the location of stars show- 
ing high rates of mass loss on the HR diagram, 
as shown schematically in Figure 5-11. Roughly 
speaking, massive stellar winds, in the range 
10' 8 to 10~ 5 M q /yr are found only among the 
more luminous stars lying in the upper half of 
the HR diagram. It is tempting to search for 
a single mechanism that can explain mass loss 
in both hot and cool stars, but this seems 
unlikely in view of the probable dependence of 
mechanisms on the environment. Thus, in the 
hot stars, the radiative force on atoms is op- 
timized by the occurrence in the same 
wavelength region of both a strong radiation 
field and the resonance lines of abundant ions, 
whereas conditions in the low-temperature stars 


are favorable for mass ejection by radiation 
pressure on dust grains. It also does not seem 
probable that the same mechanism can produce 
terminal velocities of 10 km/s in the cool stars 
and 3000 km/s in the hot stars. This section will 
be devoted to mechanisms that seem appro- 
priate for the late-type giants, and the reader 
is referred to the monographs on B stars and 
on O, Of, and Wolf-Rayet stars for an account 
of hot-star mechanisms. 

Figure 5-12, from a compilation by Reimers 
(1977b), shows the regions in the HR diagram 
in which late-type stars are found to have 
massive winds. All G and K supergiants, the 
later K giants, and all M giants and supergiants 
show circumstellar spectral lines. The region of 
the diagram containing stars with circumstellar 
dust shells is also shown, as are the locations 
of maser stars and long-period variables. IUE 



Figure 5-11. HR diagram in which dashed lines outline regions in which stars with high mass loss 
are found. Also shown are the zero-age main sequence (ZAMS), the giant branch, the locations of 
masers and long-period variables, evolutionary tracks, and the region occupied by stars showing the 
10-g.m dust signature (Goldberg, 1979). 
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Figure 5-12. Portion of HR diagram for late-type stars , showing regions in which 
stars have expanding circumstellar shells. Approximate expansion velocities are 
shown f as well as the transition zone in which winds may be variable (Cassinelli, 
1979). 


and HEAO observations reveal that all late- 
type stars have chromospheres, but only stars 
with low mass loss have transition regions and 
high-temperature coronae (Linsky and Haisch, 
1979; Vaiana et al., 1981). Linsky and Haisch 
(1979) first used IUE spectra to identify a 
boundary line in the HR diagram separating 
stars with coronae from those with only 
chromospheres. 

Stars within the zone sometimes show both 
high-temperature lines such as C IV 1550 A and 
CS lines, and the latter tend to show variabil- 
ity as well. 

Similar boundaries, although not always co- 
incident with one another, separate stars in 
which chromospheric expansion is indicated by 
the asymmetry of both the Ca II H and K emis- 
sion (Stencel, 1978) and the Mg II h and k emis- 
sion (Stencel and Mullan, 1980) and by stars 
with C IV line emission (Dupree, 1981) and X- 
ray emission (Vaiana et al., 1981). In reality, 
the boundary between stars with low and high 
mass loss is not a sharp line, but a relatively 
broad zone (Hartmann et al., 1980, 1981). 


Some stars in the boundary zone are “hy- 
brids,” so-called because their spectra show 
both displaced circumstellar lines and C IV line 
emission (Hartmann et al., 1980, 1981; 
Reimers, 1982). The CS lines in hybrid stars 
tend to show variability as well. On the average, 
terminal velocities increase from right to left, 
approaching 100 km/s near the boundary. The 
degree of sharpness of the boundary is still 
somewhat controversial (Haisch and Simon, 
1982; Simon et al., 1982). The reason for the 
boundary is undoubtedly connected with evolu- 
tionary changes creating conditions favorable 
for high mass loss, but the nature of the 
changes is still not clear. On one side of the 
boundary, nonradiative energy beyond the 
chromosphere is dissipated into heat, and on 
the other side, the energy goes mainly into driv- 
ing massive winds. 

Several types of forces have been proposed 
to account for steady mass loss in red-giant 
stars, the principal ones being: (1) thermal gas 
pressure, (2) radiation pressure on dust grains, 
(3) shock waves, primarily those arising from 
pulsations in Mira stars, and (4) Alfven waves. 
In addition, models of nonsteady or episodic 
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mass loss have been predicated on dynamical 
or pulsational instability, on convection, or on 
magnetic reconnection. Steady mass loss will 
be considered first. 


where v is the flow velocity, q is the density, 
P is the gas pressure, M* is the mass of the star, 
and 7(r) is the ratio of the radiative force to 
the gravitational force: 


Thermal Gas Pressure 

Parker’s original solar wind model has in- 
spired numerous attempts to drive winds in late- 
type giants by thermal pressure, most recently 
by Hearn (1975), Mullan (1978), and Watanabe 
(1981). Basically, these models fail because the 
thermal energy inferred from the observed low 
wind temperatures is too low to overcome 
stellar gravity (Holzer et al., 1983). The ob- 
served mass-loss rates can be achieved only by 
postulating temperatures so high that the radia- 
tion emitted would far exceed that observed. 



where Q is the flux weighted mean of the radia- 
tion pressure efficiency of the grains, n gr is the 
number density of the grains, a is the grain 
radius, and L is the stellar luminosity. 

These equations must be solved in conjunc- 
tion with the equation of motion of the grains: 


Radiation Pressure on Dust Grains 

A complete theory of dust-driven mass loss 
should include consideration of matters such 
as grain formation (see Kwok, 1975; Deguchi, 
1980; Draine, 1981; Lefevre et al., 1982; 
Woodrow and Auman, 1982), the transfer of 
infrared continuum radiation (Jones and Mer- 
rill, 1976; Menietti and Fix, 1978; Rowan- 
Robinson, 1980), circumstellar shell chemistry 
(Goldreich and Scoville, 1976; Scalo and 
Slawsky, 1980; Huggins and Glassgold, 1982; 
Glassgold and Huggins, 1983 and this volume), 
and the dynamics of the gas and dust flow 
(Kwok, 1975; Goldreich and Scoville, 1976; 
Menietti and Fix, 1978; Tielens, 1983). 

Calculations of the motion of gas in a cir- 
cumstellar shell always assume a spherically 
symmetric atmosphere with steady flow gov- 
erned by the equations of conservation of mass 
and momentum (Kwok, 1975): 
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where v gr is the velocity of the grain, q s is its 
density, and F dmg is the drag force produced 
by collisions with gas molecules. It is easy to 
show that, for grain radii of about 0.1 pm in 
late-type giants, the gravitational force is much 
smaller than the radiative force and can 
therefore be neglected. It has also been found 
that, for a Ori for example, the terminal veloci- 
ty is reached by grains in a distance on the order 
of one solar radius. Therefore, 


F = F 
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1975): 


— ip • v • r 2 ) = 0 , 
dr 


(5-15) 


F dra = a • *a 2 • 6 v d (v\ + v 2 )" 2 - (5-20) 




GM, 


[1-3 O')] 


(5-16) 


where q is the gas density, v d is the drift 
velocity of the dust grains relative to the gas 
molecules, v T is essentially the thermal veloci- 
ty, and a is an elasticity factor assumed to be 
equal to 3/4. At this point, additional equations 
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must be introduced according to whatever 
assumptions are made about the temperature 
distribution, the sites of grain formation, the 
fraction of gas in the form of grains, and their 
growth and sputtering rates. 

The assumptions that are made often de- 
pend on the primary goals of the models. Thus, 
the models of Kwan and Hill (1977), Kwan and 
Linke (1982), and Morris (1980), which were 
described in the section Millimeter Waves , 
assumed a priori that the mass loss is driven by 
dust and were used to calculate molecular line 
emission, with the rate of mass loss as an ad- 
justable parameter to be determined from ob- 
servation. In this section, we shall be concerned 
with models whose purpose it is to test various 
aspects of the dust-driven mechanism. 

Following the approximate analytical ex- 
pressions derived by Gehrz and Woolf (1971), 
Kwok (1975) made the first numerical solutions 
of the equations of motion, taking account of 
radiation pressure on grains, stellar gravity, 
growth and sputtering of grains, and momen- 
tum transfer from grain to gas. Kwok believed 
that radiative heating and cooling could be 
neglected and assumed the adiabatic law for the 
temperature variation. In their investigation of 
the physical properties of OH/IR stars, Gold- 
reich and Scoville (1976) showed that collisions 
between gas molecules and dust grains are a 
significant source of heating and that the gas 
is cooled both by adiabatic expansion and by 
H 2 0 molecular emission. The resulting tem- 
perature distribution is markedly different from 
the adiabatic law. Goldreich and Scoville also 
considered the chemical reactions that control 
the abundances of OH and H 2 0. 

Menietti and Fix (1978) solved the equation 
of motion of the gas in cool luminous stars, tak- 
ing account of the radiative properties of both 
the gas and the dust. The density, temperature, 
and radiative energy distributions were not 
assumed in advance, but were derived from the 
modeling process. Similarly, Rowan-Robinson 
(1980) has developed accurate solutions to the 
equation of radiative transfer for the purpose 
of predicting the radiative energy distributions 
of stars throughout the wavelength region 0.4 


to 30 fim . The model consists of a blackbody 
of given radius and temperature at the center 
of a spherically symmetric dust shell. The free 
parameters of the model are the temperature 
of the central star, the condensation tempera- 
ture of the grains, which is also the temperature 
of the inner boundary of the shell, the limiting 
optical depth at short wavelengths through the 
shell, the density gradient in the shell, the ratio 
of the inner and outer radii of the shell, and 
the interstellar extinction. The parameters are 
determined by fitting model calculations to the 
observed spectra. Models have been derived 
and applied to the observed spectra of 156 late- 
type stars, including early M stars (Rowan- 
Robinson and Harris, 1982), stars of type M5 
and later (Rowan-Robinson and Harris, 1983a), 
carbon stars (Rowan-Robinson and Harris, 
1983b), and OH/IR stars showing double shells 
(Rowan-Robinson, 1982). 

The model by Deguchi (1980), of the flow 
in O-rich stars, focuses mainly on the process 
of grain formation. The theory of grain nuclea- 
tion is used to calculate the number density of 
grains, the grain sizes, and the final amount of 
metals remaining in gaseous form. Woodrow 
and Auman (1982) construct time-dependent 
models of carbon-rich red giants, making use 
of the equations of hydrodynamics, the theory 
of grain nucleation and growth, and the equa- 
tions describing a gray model atmosphere in 
LTE and radiative equilibrium. 

The most recent model of dust-driven gas 
flow in O-rich Mira stars (Tielens, 1983) is also 
the most detailed and complete. The starting 
point of the model is the observational evidence 
(Hinkle, et al., 1982) that the outward motion 
of the pulsating photospheres of Miras gener- 
ates shock waves that elevate matter to a sta- 
tionary layer about 5 stellar radii above the 
photosphere, where the temperature is less than 
1000 degrees and therefore low enough to per- 
mit dust grains to condense. The stationary 
layer is, in effect, a reservoir for the matter in 
the wind. 

Some general conclusions can be drawn 
from these modeling investigations. Cir- 
cumstellar gas shells in late-type giants and 
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supergiants are heated by stellar radiation and 
by grain/gas collisions. At least in the late-M 
Miras and supergiants, the gas is heated by the 
drag force and cooled by H 2 0 emission and 
adiabatic expansion (Goldreich and Scoville, 
1976). The velocity gradient also plays an im- 
portant part in the determination of the 
temperature structure (Tielens, 1983). When the 
density is high enough, photon trapping can 
slow radiative cooling in the absence of a veloc- 
ity gradient. Thus, the temperature structure of 
the gas is coupled to the dynamics of the out- 
flow. 

Observed radiative energy distributions are 
consistent with an r ~ 2 law for the density 
distribution (Rowan-Robinson and Harris, 
1982, 1983a, 1983b). However, detailed model- 
ing of the flow (Menietti and Fix, 1978; Tielens, 
1983) shows that the density falls off much 
more rapidly than as r ~ 2 until after the sonic 
point is reached. Correspondingly, in Tielens’ 
model, the velocities increase rapidly from the 
stationary layer, at distance r., to a terminal 
value at about 10 r. An example given by 
Tielens (1983) shows that, if the mass-loss rate 
is calculated from the total mass column den- 
sity by Equation (5-3) on the assumption that 
the velocity is constant, the rate will be over- 
estimated by a factor of 3 to 4. 

The model calculations also demonstrate 
that it is difficult to produce 10-jun excesses in 
O-rich stars with clean silicates (Jones and Mer- 
rill, 1976; Menietti and Fix, 1978; Rowan- 
Robinson and Harris, 1982, 1983a). In the car- 
bon stars, amorphous carbon grains with ab- 
sorption efficiency proportional to frequency 
give better agreement between observed and 
calculated energy distributions than does 
graphite (Rowan-Robinson and Harris, 1983b). 
Empirical values of the optical depth of the 
silicate feature at 10 /im, derived from com- 
parison of model calculations with observation, 
are strongly model-dependent. For example, 
the values derived by Rowan-Robinson and 
Harris (1982, 1983a) are smaller than those of 
Hagen (1978), in some cases by a factor of 20 
and in others by a factor of about 200. 


Differences in assumed grain parameters can 
account for a factor of 20, whereas another fac- 
tor of 10 may occur in stars for which the use 
of an optically thin model (Hagen, 1978) may 
not have been justified. Until these disagree- 
ments can be resolved, it is premature to draw 
conclusions from dust column densities inferred 
from the 10-/un optical depth (Hagen, 1978; 
Hagen et al., 1983; Bernat, 1981). 

The fraction of metals remaining in the 
gaseous state varies with the mass-loss rate 
(Deguchi, 1980) from 1 to 10 percent for M = 
10“ 5 M o /yr, and < 1 percent for 3 x 10' 6 
Af © / yr. Owing to the large drift velocity of the 
grains relative to the gas, the dust-to-gas ratio 
is not a constant with distance (Tielens, 1983). 
This is important in the estimation of mass-loss 
rates from the 10-^m feature. 

We come finally to the most important ques- 
tion. Does radiation pressure on dust grains 
drive the observed mass flows in late-type giants 
and supergiants? The most important param- 
eters on which the success of the mechanism 
depends (Kwok, 1975; Menietti and Fix, 1978) 
are the stellar luminosity, L, the mass, M*, 
which together determine the distance of the 
sonic point, the effective temperature, T ff , 
and the most critical factor, the density of gas 
at the condensation point, which must be 
large enough both to ensure momentum cou- 
pling and to hold down the sputtering rate by 
keeping the drift velocity low. The condensa- 
tion-point gas density is a function of both ef- 
fective temperature, T ff , which controls the 
distance of the condensation point from the 
star, and the scale height, which depends on the 
stellar mass, or luminosity, L. Kwok (1975) has 
calculated the minimum density at the conden- 
sation point required for mass ejection from 
stars of various luminosities and effective 
temperatures, and has shown that the more 
luminous the star, the lower the temperature at 
which mass loss is possible. Figure 5-13 is a plot 
of log L versus T max , the maximum effective 
temperature for which dust-driven mass loss 
can occur. Also plotted are the positions in the 
diagram of four red supergiants, a Ori, a Sco, 
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Figure 5-13. The maximum effective tem- 
perature at which dust-driven mass loss would 
be expected to occur , as a function of stellar 
luminosity according to model calculations by 
Kwok (1975). The positions of the four super- 
giants, a Ori , a Sco, a Her , and y Cep , are in- 
dicated by crosses. 

a Her, and fi Cep. Only for a Her is the tem- 
perature low enough to allow condensation at 
high enough gas densities. These calculations 
are not necessarily fatal to the dust-driven 
hypothesis for the hotter stars (Draine, 1981). 
For example, the gas density may be increased 
by clumping or by the ejection of discrete 
clouds, and “clean” grains can form relative- 
ly close to the star and gather up impurities as 
they move away. Note also that in a Ori, and 
probably in other red giants as well, the 
photosphere is greatly distended, its scale height 
being on the order of a stellar radius (Balega 
et ah, 1982). Calculations based on plane- 
parallel model atmospheres give erroneously 
high temperatures for spherically extended at- 
mospheres (Schmid-Burgk and Scholz, 1981). 
It remains to be investigated whether the 
temperatures are lowered sufficiently in the 
higher density regions to allow grain formation. 

An empirical test of the dust mechanism can 
be made by comparing observed rates of mass 
loss with rates estimated on the assumption that 
the total energy radiated by grains in the in- 
frared equals the energy absorbed in imparting 


momentum to the grains (Knapp et al., 1982). 
The rate of grain-driven mass loss is related to 
the infrared flux at the Earth by 

M V 0 /D 2 = 0 F/c , (5-21) 

where /? is a factor measuring the amount of 
dust in the envelope, D is the distance to the 
star, V Q is the terminal velocity, and F is the 
bolometric flux measured at the Earth, derived 
by integrating the continuum spectrum inferred 
from published IR fluxes. On the premise that 
the dust mechanism works for the oxygen stars, 
Knapp et al. (1982) calculate (3 from Equation 
(5-21) and find it to be correlated with the color 
(3.5 fim) - (1 1 yi n), which is known to measure 
the dust content of circumstellar envelopes. The 
correlation may be the best evidence to date 
that the dust mechanism is responsible for mass 
loss in the cooler O-rich stars. It is still not clear 
whether or not the dust mechanism is adequate 
to drive the mass loss in the carbon stars or 
whether some additional mechanism needs to 
be invoked. 

Shock Waves 

The conclusion from both theoretical and 
observational considerations is that radiation 
pressure on dust grains can drive mass loss at 
the observed rates in late M and S giants and 
supergiants, but that some additional force or 
forces are needed in the early M-types and 
especially in the carbon stars. In addition to 
radiation pressure on grains and thermal pres- 
sure, wave pressure has received considerable 
attention as an important source of momentum 
and energy in late-type stars, either from shock 
waves generated by pulsation in Mira stars 
(Willson and Hill, 1979; Willson, 1981; Wood, 
1979, 1981) or from Alfven waves (Hartmann 
and McGregor, 1980, 1982; Hartmann and 
Averett, 1984). Calculations of the dynamic 
response of model Mira atmospheres to a series 
of periodic shocks driven by an oscillating 
piston have shown that high rates of mass loss 
are possible in such models. Two extreme 
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cases have been studied. Both Willson and Hill 
(1979) and Wood (1979) found that isothermal 
shock models showed very low rates of mass 
loss, about 10" 12 M q /yr, whereas in the 
adiabatic limit, the shocks give unrealistically 
high rates, about 10 2 M© /yr, and unaccept- 
ably large values of the terminal velocity. Thus, 
the critical point in the theory is whether the 
shocks generated by pulsation are isothermal 
or adiabatic. Some compromises are possible: 

1 . Periodic isothermal shocks can levitate 
the atmosphere and thereby increase 
grain-driven mass loss by a factor of 
about 40, according to Wood (1979). 
Another estimate by Jones et al. (1981) 
is that the increase in atmospheric scale 
height by pulsation can enhance the rate 
of grain-driven mass loss by several 
orders of magnitude. 

2. In the lower atmosphere, where the 
emission lines originate, the shocks may 
be isothermal, but higher up, as the den- 
sity falls, the rates of cooling diminish 
and the shocks may become adiabatic 
(Willson and Hill, 1979). 

Alfven Waves 

Hartmann and McGregor (1980, 1982) have 
studied the response of model atmospheres of 
red giants and supergiants to the passage of 
Alfven waves of about solar magnitude, with 
stellar gravity as the main parameter. They con- 
clude that the force due to Alfven waves can 
account for observed mass-loss rates and ter- 
minal velocities, provided that the waves are 
dissipated over a distance of the order of a 
stellar radius by damping arising from collisions 
between ions and neutral atoms. The theory is 
an attractive one, both because a single 
parameter — stellar gravity — controls the rate of 
mass loss and because it predicts that red giants 
have extended warm chromospheres as ob- 
served. However, a detailed study of Alfve'n 
waves in stellar atmospheres (Holzer et al., 


1983) concludes that cool winds driven by 
Alfven waves will not exhibit both high mass 
loss and low terminal velocity. The basic prob- 
lem with this and other stellar wind mech- 
anisms, as Holzer et al. (1983) point out, is that 
energy inserted into the flow in the subsonic re- 
gion increases the mass flux but not the veloc- 
ity, whereas in the supersonic region, the addi- 
tion of energy increases the terminal speed 
without affecting the mass loss. Thus, the 
mechanism must be such that a large fraction 
of the energy is dissipated in the subsonic flow 
region, with just the right small amount left 
over for the supersonic flow to overcome gravi- 
ty and yet limit the terminal velocity to the 
observed fraction of the velocity of escape from 
the photosphere. 

Recently, Hartmann and Averett (1984) 
have applied the theory to construct a model 
of the extended chromosphere of a Ori. Assum- 
ing that a wind of about 10" 6 M© /yr is driven 
by dissipating Alfven waves, they reproduce the 
observed free-free emission at radio wave- 
lengths, the Ha emission observed by speckle 
interferometry, and the emission fluxes in a 
number of spectral lines. However, the model 
fails to reproduce the observed shapes of 
chromospheric line profiles, which suggests 
that, although the waves can heat and extend 
the chromosphere, the observed chromospheric 
expansion velocities are three times smaller than 
those predicted by the model. Hartmann and 
Averett (1984) suggest two possible ways out 
of this problem. Either the waves simply dis- 
tend the atmosphere so that the gas density at 
the grain condensation point is high enough to 
maintain dust-driven mass loss, or spatial and 
temporal time variability in wave generation 
and damping may lead to composite emission 
from low- and high-velocity ejections. 

If massive winds in late-type stars are driven 
by Alfven waves, the stars should have extend- 
ed chromospheres, which emit microwave radi- 
ation at centimeter wavelengths. Such radi- 
ation has been detected from several K and 
early M giants and supergiants (Drake and Lin- 
sky, 1984), but its weakness or absence in 31 


278 


evolved red giants surveyed by Spergel et al. 
(1983) argues against the effectiveness of the 
Alfven-wave mechanism in the late M stars. 

Mechanisms for Nonsteady Mass Loss 

A number of mechanisms have been pro- 
posed for sporadic or nonsteady mass loss in 
red giants. Schwarzschild (1975) suggested that 
large-scale convection might trigger the ejection 
of a massive cloud by some process as yet un- 
known. Such ejections might cause both the ir- 
regular brightness variations observed in cool 
variable stars and the observed ordered changes 
in the linear polarization. Bloemhof et al. 
(1984) suggest that such localized and tempo- 
rary phenomena are consistent with the ob- 
served asymmetric distribution of dust around 
a Ori. Wood (1979) has proposed that isother- 
mal shock waves in the atmosphere of a Mira 
star generate mass loss by levitating the atmo- 
sphere to heights at which grains condense. 
Jones et al. (1981) have shown that the mass- 
loss rate may be enhanced by several orders of 
magnitude, as compared with grain-driven mass 
loss alone. We have already referred to the evi- 
dence that such a process may be at work in 
X Cygni and in other Miras as well (Hinkle et 
al., 1982). In stars like a Ori, matter might be 
transported to the dust layer by shock waves 
associated with the irregular fluctuations in 
brightness, as suggested, for example, by 
Schwarzschild (1975) and by Kafatos and 
Michalitsianos (1979). 

The region in a Ori between the photosphere 
and the base of the shell, which we tentatively 
identify as the chromosphere, may be the 
source of the Fe II emission lines observed by 
Boesgaard and Magnan (1975) and Boesgaard 
(1979) near 3100 A, by van der Hucht et al. 
(1979) in the 2750 to 3165 A region, and by Car- 
penter (1984) in the 2300 to 3000 A region. 
Measurements of these lines have given con- 
flicting results. Boesgaard and Magnan found 
that the motions derived from the Fe II lines 
were closely coupled to the photospheric mo- 
tions, but indicated matter falling into the 
photosphere at an average speed of about 


5 km/s. Conversely, van der Hucht et al.’s mea- 
surements from a balloon showed an average 
blue-shift of 14 km/s relative to the photo- 
sphere. From high-resolution IUE spectra, 
Carpenter (1984) found no evidence of strong 
coupling between photosphere and chromo- 
sphere and no systematic difference between the 
mean velocities of the photosphere and chro- 
mosphere. He interpreted the relation between 
line asymmetry and line strength as evidence for 
a schematic velocity field that is initially con- 
stant with height and then increases outward 
to a maximum value before decreasing to an 
asymptotic value. These apparently contradic- 
tory results might be resolved if the velocity 
fields varied with time, which might be expected 
if the Fe II emission is responsive to fluctua- 
tions in photospheric activity. On the other 
hand, Carpenter’s data are more complete and 
of higher quality than those obtained from the 
ground and with balloons. Continued monitor- 
ing from IUE would be very important in re- 
solving this question. 

Burke (1972) suggested that rising shock 
waves ejecting matter would be expected when- 
ever there is pulsation or a strong convection 
zone near the surface of a star. The mass would 
be ejected preferentially from the equator in 
the direction of rotation and would lead to loss 
of angular momentum without the need for 
magnetic fields. Kafatos and Michalitsianos 
(1979) have applied this mechanism to the de- 
termination of the rate of spin-down of red su- 
pergiants. The ejection of matter is predicted 
to form a disk or torus of silicate grains in the 
equatorial regions. 

Mass loss is also predicted to occur from dy- 
namically and pulsationally unstable envelopes 
of red-giant stars (Tuchman et al., 1978, 1979, 
and earlier references therein). According to 
Tuchman et al. (1978), dynamical instability 
sets in when, owing to partial ionization and/or 
partial dissociation, the quantity, 7 , the 
adiabatic value of the logarithmic pressure gra- 
dient, falls below 4/3 over a sufficiently large 
and/or deep region. An initial expansion turns 
into violent oscillations in the envelope. Shocks 
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generated within the envelope lead to repetitive 
episodes of mass loss on a time scale of < 30 
years, which would cause large amounts of 
mass to be lost during evolutionary time scales. 
Newer calculations by Fox and Wood (1982) 
show that the apparent instabilities are caused 
by the neglect of nonadiabatic effects in the 
models. 

Magnetic Reconnection 

A few years ago, Mullan (1978) suggested 
that the boundary in the HR diagram between 
stars with high and low mass loss is a super- 
sonic transition locus (STL), along which the 
sonic points of the wind dip below the heights 
reached by the tops of spicules. Holzer (1980) 
and others have shown the STL hypothesis to 
be untenable, and Mullan (1982) has proposed 
instead that the boundary is caused by a trans- 
ition in magnetic topology. (See also Mullan 
and Ahmad, 1982; Mullan and Steinolfson, 
1983). The idea is that, since X-ray emission in 
the Sun is associated with closed magnetic loops 
and the solar wind is enhanced over coronal 
holes where the field is mainly open, massive 
winds and the absence of X rays may result 
when stars are unable to maintain closed mag- 
netic field lines in the steady state. With the aid 
of Pneuman’s (1968) model of helmet stream- 
ers, Mullan (1982) investigates the question of 
where in the HR diagram the transition from 
closed to open field configurations might oc- 
cur. Pneuman (1968) found that helmet stream- 
ers would be unstable and would not exist in 
a steady state when the parameter, 

* = n - G • A/, • m p /R,kT , (5-22) 

is smaller than 4. In Equation (5-22), /x is the 
mean molecular weight, m p is the proton 
mass, M, and are the mass and radius of 
the star, respectively, and T is the temperature 
at the sonic point. Since the radius, R s , of the 
sonic point is \lGMjyiJ 2 kT , the condition 
that long-lived closed lines can exist is R $ > 2 
R m . The locus R s = 2 R* on the HR diagram 
defines a so-called magnetic topology transition 


locus (MTTL) which happens to closely follow 
the boundary between stars with weak and 
strong mass loss. Once it is accepted that the 
onset of massive winds in red giants is 
associated with open field lines, the Alfven- 
wave mechanism may be invoked to drive the 
wind. Alternatively, Mullan (1981) has spec- 
ulated that after unstable emerging loops are 
ejected from the surface, field-line reconnec- 
tion at the base of the loops furnishes the 
energy for accelerating them upward. 

EVOLUTIONARY CONSEQUENCES 
OF MASS LOSS 

Supernovae 

The evolutionary effects of mass loss on red 
giants have been well reviewed, for example, 
by Iben (1981b) and Renzini (1981a, 1981b). 
The most convincing evidence comes from 
supernovae. In the absence of mass loss, all 
stars with initial masses less than 1.44 M© 
would become white dwarfs, whereas all those 
with larger masses would explode as supernovae 
of either Type I or Type II, depending on the 
mass. It is estimated that, without mass loss, 
the frequency of supernovae in a galaxy would 
be one per year, which is about 10 times greater 
than observed. Observed rates of mass loss 
allow one to predict that stars with initial 
masses as large as 5 to 8 M© can end up as 
white dwarfs. The predicted frequency of 
supernovae is reduced to about the observed 
level. 

Stars in Globular Clusters 

A second theoretical argument for mass loss 
in red giants is morphological and is based on 
the observed distribution of stars on the giant 
and horizontal branches of Population II 
globular clusters. Evolutionary calculations 
without mass loss give perfect fits to the main- 
sequence and red-giant branches of the HR 
diagram, but the tip of the asymptotic giant 
branch (AGB) is about 2.5 mag too bright and 
the horizontal branch (HB) is too short: blue 
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stars and RR Lyrae stars are absent. Stars em- 
barking on the red-giant branch (RGB) are 
estimated to have initial masses of 0.8 to 0.9 
M o . The HB discrepancy may be removed if 
stars lose 0.2 M © solar masses during the RGB 
phase before helium flash, whereas the AGB 
problem is resolved if an additional 0.1 M© 
is lost during the AGB phase. The computed 
HB tracks are still too short, unless it is as- 
sumed that stars arriving at the HB have a 
spread in their masses of about 10 percent. 

Evolutionary calculations incorporating 
mass loss as given by the revised Reimers form- 
ula, Equation (5-4), reproduce the observed 
morphology of the RGB, HB, and AGB in 
Population II globular clusters. It is claimed 
further that such morphological considerations 
put narrow constraints of about 10 percent on 
the correct values of mass-loss rates themselves. 
Judgment on this point must be reserved until 
more is known about the interior physics of 
stars (e.g., convection, neutrino physics, tur- 
bulent diffusion, primordial chemical composi- 
tion, and meridional mixing). 

Composition of Interstellar Matter 

The theory of stellar evolution also predicts 
that red supergiants may be the most impor- 
tant source of interstellar grains, both from 
mass loss and from supernova explosions. 
Variations in the surface composition during 
the AGB phase should affect grain character- 
istics in the winds. At the beginning of the AGB 
phase, the star is type M with C/O < 1 , and con- 
sequently, silicate grains can form in the wind. 
Later, the star becomes a carbon star with C/O 
>1; after the third dredge-up (Iben, 1981a), O 
is locked up in CO, but there is an abundance 
of free carbon and C-rich molecules. Hence, 
grains in the winds will contain not silicates but 
ample quantities of C-rich particles such as 
amorphous carbon grains, graphite, and SiC 
grains. 

Formation of Planetary Nebulae 

The formation of planetary nebulae is one 
of the most spectacular consequences of mass 


loss from red giants, and in fact, observations 
of planetary nebulae are providing critical tests 
of the entire theory of the late stages of stellar 
evolution, in which mass loss plays a critical 
role. It is generally accepted that planetary 
nebulae are produced by asymptotic branch 
stars, which eject their hydrogen-rich enve- 
lopes, after which the stellar remnant evolves 
to a high temperature and radiates the far-UV 
photons needed to excite the envelope. By what 
mass-loss process is the envelope ejected and 
at what rate? Assuming with Renzini (1981b) 
that a typical planetary nebula has a mass of 
0.2 Mq , a radius of 0.1 pc, and an expansion 
velocity of 20 km/s, we find that the nebula 
must form in about 5000 years, which requires 
a mass-loss rate of about 4 x 10' 5 Mq/ yr. 
Since this rate is an order of magnitude greater 
than that given by the Reimers formula when 
applied near the top of the AGB, Renzini 
postulated the more or less sudden onset of a 
“superwind” that would turn on at the end of 
the AGB phase and operate for a relatively 
short time. 


Such a superwind could take one of several 
forms. It could occur when Mira variables 
switch oscillation modes from the first overtone 
to the fundamental (Wood, 1974, 1981; Tuch- 
man et al., 1979; Fox and Wood, 1982) or by 
the onset of pulsation in Miras, which Willson 
(1981) believes is normally in the fundamental 
mode. It was once believed that the formation 
of planetary nebulae might be caused by dy- 
namical instability, but detailed calculations 
(Tuchman et al., 1978) have shown that the 
relevant models should stabilize while the 
envelope mass is still too large (> 0. 1 Mq ). The 
mode-switching mechanism leads to a relative- 
ly sudden ejection of the stellar envelope in less 
than 1000 years (Tuchman et al., 1979). A dif- 
ferent scenario which does not require sudden 
ejection has been proposed by Kwok et al. 
(1978) and has been elaborated by Kwok (1982). 
The steady mass loss eventually depletes the H- 
rich envelope, exposing the hot core and caus- 
ing a sudden change in the mass-loss mech- 
anism to one of radiation pressure on ions. The 
velocity of the wind increases by a factor of 
100, and the snowplow effect builds up a dense 
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shell at the interface between the new and the 
old wind. In this model, the planetary nebula 
forms in about 3000 years. The distinction be- 
tween steady and sudden ejection seems arti- 
ficial, inasmuch as both mechanisms require a 
sudden increase in mass-loss rate, and the dif- 
ference in time scales is well within the uncer- 
tainties of the calculations. Moreover, interac- 
tion is bound to occur between the old circum- 
stellar shell and the new wind, no matter which 
mechanism is operative. 

Recent observations of unidentified OH/IR 
stars (Baud et al., 1981; Engels et al., 1983) are 
providing considerable evidence in favor of 
Mira variables as progenitors of planetary 
nebulae, as was suggested by Wood and Cahn 
(1977). As stars evolve up the AGB, they lose 
mass steadily, and at some value of the lumi- 
nosity, which depends on the mass, they begin 
pulsating as Mira variables in the first-overtone 
mode. As the luminosity increases, so do the 
radius and the period. Eventually, for a given 
stellar mass, a star reaches a value of the 
luminosity at which the fundamental mode be- 
comes dominant, the envelope is ejected, and 
the planetary nebula forms in a matter of a few 
thousand years. The evolutionary model of 
Wood and Cahn (1977) was developed in con- 
siderable detail by Tuchman et al. (1979), but 
encountered several disagreements with obser- 
vation, the most serious being the termination 
of first-overtone pulsation at much lower lu- 
minosities and shorter periods than observed, 
and the occurrence of dynamical instabilities at 
the point at which the fundamental mode takes 
over, which limits the observed periods to less 
than 600 days. It now appears that the insta- 
bilities do not in fact occur (Fox and Wood, 
1982) and hence that pulsations with periods 
longer than 2000 days are stable. 

From a study of the statistical properties of 
about 20 OH/IR stars, Baud and Habing (1983) 
were able to investigate the physical parameters 
that determine the OH luminosity. The OH 
luminosity is correlated with the size of the cir- 
cumstellar shell and with the mass-loss rate. 
They make the basic assumption that the se- 


quence of OH/IR stars arranged in order of in- 
creasing period and mass-loss rate (see Table 
5-3) is an evolutionary sequence, and they 
derive a timescale of between 0.6 x 10 5 and 
6.5 x 10 5 years for main-sequence masses be- 
tween 0.7 to 6 M © for the OH maser lifetime. 
Of this time, 99 percent is spent as a visible 
Mira with low mass loss, and 1 percent is spent 
as a luminous OH/IR star. Comparison of the 
derived physical parameters of the stars with 
the models of Fox and Wood (1982) suggests 
that the OH maser turns on at about the time 
that the pulsation mode switches from first 
overtone to fundamental. The models also pre- 
dict that the longest periods go with the highest 
mass-loss rates, as found by Engels (1982). 
Engels et al. (1983) also conclude that OH/IR 
stars are AGB variables, but that they are prob- 
ably in the same evolutionary stage as Miras. 
They argue that, since OH/IR stars show higher 
concentration to the galactic plane than Miras, 
they must have higher masses, in the range 2 
to 10 Mq. The high galactic concentration of 
unidentified OH/IR stars is not a problem 
(Baud and Habing, 1983), inasmuch as these 
stars represent only 1 percent of the OH/IR and 
are distant objects. It is less clear why there 
should be no low-mass OH/IR stars with high 
OH luminosity. 

Radio continuum emission at 6 cm, which 
may be emitted by ionization fronts propa- 
gating into CS shells during the preplanetary 
nebula stage, has been detected in five of 32 
cool evolved stars: NGC 7027, CRL 618, R 
Aqr, IRC + 10216, and O Ceti (Spergel et al., 
1983). Model calculations predict that, in ob- 
jects with low mass-loss rates, M < 10“ 
A/ 0 /yr, a planetary nebula is produced in 
< 10 3 years after the hot core is exposed. In 
those with intermediate masses, 4 to 6 M Q , and 
high mass-loss rates, the ionization front moves 
so slowly that the shell is always ionization- 
bounded, and most of the mass is ejected as 
neutral gas and dust. NGC 7027 is an example 
of a recently formed planetary nebula, while 
CRL 618 is a protoplanetary nebula with its hot 
core still obscured by dust, and IRC + 10216 
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may or may not ever reach the planetary stage. 
Thus, the requirement for PN formation is a 
short period of very high mass loss, such as is 
found in long-period OH/IR sources, followed 
by one of reduced mass loss to allow the rapid 
propagation of an ionization front. 
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